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Abstract 

This paper describes perturbative framework, on the basis of closed-time-path formalism, for studying quasiuniform 
relativistic quantum field systems near equilibrium and nonequilibrium quasistationary systems. At the first part, 
starting from first principles, we construct perturbative schemes for relativistic complex-scalar-field theory. We clarify 
what assumption is involved in arriving at a standard perturbative framework and to what extent the n (> 4)-point 
initial correlation functions that are usually discarded in the standard framework can in fact be discarded. Two 
calculational schemes are introduced, the one is formulated on the basis of the initial particle distribution function 
and the one is formulated on the basis of the "physical" particle distribution function. Both schemes are equivalent 
and lead to a generalized relativistic kinetic or Boltzmann equation. At the second part, using the perturbative loop- 
expansion scheme for an O(N) linear a model, we analyze how the chiral phase transition proceeds through disoriented 
chiral condensates. The system of coupled equations that governs the spacetime evolution of the condensate or order- 
parameter fields is derived. The region where the curvature of the "potential" is negative is dealt with by introducing 
the random-force fields. Application to simple situations is made. 
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I. INTRODUCTION 



Since mid-fifties, efforts have been made to incorporate quantum field theory with nonequilibrium statistical me- 
chanics ji]-f|]. Necessity of this incorporation originated from the field of solid-state physics. Since then, rapid 
progress of the studies of the early Universe and the quark-gluon plasma have further activated this field of research 
(cf., e.g., [[||(|). The quark-gluon plasma is a system of (anti)quarks and gluons, which is governed by quantum 
chromodynamics (QCD). According to QCD, there are two types of phase transitions; a chiral phase transition and a 
confined-deconfined phase transition. According to lattice analyses, critical points for these two transitions are very 
close. The quark-gluon plasma is the system that is in the deconfined and chiral-symmetric phases and is expected to 
be produced in relativistic heavy-ion collisions and to have existed in the early Universe. Such quark-gluon plasmas 
evolve into a confined (or hadronic) and chiral-symmetry breaking phases through the phase transitions. The study of 
this evolution on the basis of QCD is very involved and some tractable effective theories have been proposed. Among 
those is a linear a model, which was proposed long ago ]tJ and is revived as a toy model that mimic the aspect of the 
chiral phase transition of QCD]]. A numerous work was reported so far on the "static" properties of this model or its 
variants. (See, e.g., ]9|Jl0|. Earlier work is quoted therein.) Recently, analyses of dynamical evolution of the system 
from a symmetric phase to a broken symmetric phase have begun (see, e.g., pd[). Such an analysis is still involved 
and different approximations have been employed by various authors. 

The purpose of this paper is twofold. In the first part, we aim at laying perturbative framework for dealing with 
out-of-equilibrium relativistic quantum complex-scalar-field systems. In the second part, taking up an fermionless 
O(N) linear a model, we lay down the framework, based on the perturbative loop expansion of the effective action, 
for analyzing how the phase transition proceeds through disoriented chiral condensates (see, e-g;, [|l2]]). We use 
the standard framework of nonequilibrium statistical quantum-field theory that is formulated — |3|,pL3|| by employing 
the closed-time path, — oo — > +00 — > —00, in a complex-time plane, which is referred to as the closed-time-path 
(CTP) formalism. Specializing the density matrix, which characterizes the initial states of the systems, to the one 
for the Gibbs distribution, the framework reduces to the standard framework of equilibrium thermal quantum field 
theories [^,0. Throughout this paper, we are interested in quasiuniform systems near equilibrium or nonequilibrium 
quasistationary systems, which we simply refer to as out-of-equilibrium systems. This means, in particular, that the 
phase transition dealt with in the second part proceeds slowly. 

The out-of equilibrium systems are characterized by two different spacetime scales: microscopic or quantum-field- 
theoretical and macroscopic or statistical. The first scale, the microscopic-correlation scale, characterizes the range of 
radiative correction to reactions taking place in the system while the second scale measures the range of interaction of 
particles. For a weak-coupling theory, in which we are interested in this paper, the former scale is much smaller than 
the latter scale. A well-known intuitive picture (cf., e.g., fl3|| ) for dealing with such systems is to separate spacetime 
into many "cells" whose characteristic size, L M (fi = 0, ...,3), is in between the microscopic and macroscopic scales. 
It is assumed that the correlation between different cells is negligible in the sense that microscopic or elementary 
reactions can be regarded as taking place in a single cell. On the other hand, in a single cell, relaxation phenomena 
are negligible. Thus, in dealing with a reaction, all the vertices in a Feynman diagram for the reaction are in a single 
cell. In more precise, propagators are Fourier transformed over a cell and carry an "index" that labels the cell. 

Above intuitive picture may be implemented as follows. Let A(x, y) be a generic propagator or two-point function. 
An inverse Fourier transformation with respect to the relative coordinates x — y yields 

A{X; P)= Jd 4 (x~ y) e «M*-iO A(x, y), (1.1) 

where X — (x + y)/2 is the midpoint, which is chosen as a label fo r th e cells. Sometimes, X is called the macroscopic 



spacetime coordinates. It is obvious from the above argument that (1.1) is meaningful for |P M | > 1/LA The self-energy 
part T,(x,y) enjoys similar property. 

Development of the phase transition is described by an evolution equation for a set of condensate or order-parameter 
fields 0{x). As mentioned above, we assume that the phase transition proceeds slowly. By this we mean that the 
x dependence of (p{x) is weak, or more precisely, x of (p(x) may be treated as a macroscopic spacetime coordinates. 
It should be noted here that, as the system approaches the critical point of the phase transition, the microscopic 
correlation scale diverges. This means that the formalism applies to the systems away from the critical region. 
Recalling that the phase transition occurs when the squared effective mass, M 2 (x), of the quasiparticles becomes 



1 Lattice simulations of finite-temperature QCD indicates that the chiral phase transition is of second order (see, e.g.. 
The phase transition described by the linear a model is of second order. 
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zero, we can say that the formalism may be used in the region, where |M 2 (cc)| is not small as compared to the typical 
scale parameter (s) of the system. 

Microscopic reactions discussed above cause change in the density matrix, through which the number densities of 
quasiparticles change with macroscopic spacetime X^. Dealing with this is the subject of the "next stage," where 
(weak) ^''-dependence of A(X;P), £(X;P), and tp(X) are explicitly taken into account. 

As mentioned above, much work has already been devoted to the issues settled above. However, most work devoted 
to the first issue (to be dealt with in the first part of this paper) assume (explicitly or implicitly) special forms for 
the density matrix. Also as mentioned above, the second issue has so far been studied under various approximations. 
In this paper, starting from first principles, we proceed with discussing matters in a context which is as general as 
possible. 

The plan of the paper is as follows: In Section II, a brief introduction is given to the CTP formalism for out- 
of-cquilibrium quantum field theories. In Sections III - V, taking up a self- interacting relativistic complex-scalar 
quantum field theory, perturbative calculational schemes are constructed. In Section III, standard representation 
of the CTP Green functions in the interaction picture is derived from first principles. In the course of derivation, 
we clarify what assumption is necessary in arriving at the representation. In a standard perturbative scheme, the 
n (> 4)-point initial correlation functions are also usually discarded. We discuss to what extent they can really be 
discarded. Whenever necessary, they can be incorporated into the scheme (cf. 0). In Section IV, we deal with the 
bare propagator (two-point function), which is a building block of perturbation theory. The quasiparticle fields are 
introduced through a sort of Bogoliubov transformation, which is a generalization of the so-called thermal Bogoliubov 
transformation in equilibrium thermal field theory |^J^JT^|. Quasiparticle fields describe, in a sense, stable modes in 
the system. With the aid of the quasiparticle picture, we obtain a 2 x 2 matrix representation for the bare propagator. 
One of the advantages of this picture is that the elementary-reaction part and the relaxation part, mentioned above, 
are clearly identified. At the final part of Section IV, we propose two perturbative schemes, which are equivalent 
to each other. In Section V, through analyzing the structure of 2 x 2 self-energy matrix, we find the form for the 
self-energy-part resummed or full propagator. Then, after discussing how the two perturbative schemes work, we 
show that both schemes are equivalent and lead to a generalized relativistic kinetic or Boltzmann equation. Sections 
VI - VIII are devoted to the analysis of a fermionless O(N) linear a model, with the aid of the effective action. 
[The case of N = 4 is of practical interest.] Here, the condensate or order-parameter fields and the quantum fields 
come in. The latter describe quantum fluctuations around the former. The condensate fields, which define the 
internal reference frames for the quantum fields, change in spacetime. Taking this fact into account, in Section VI, we 
construct consistent perturbative loop-expansion scheme for studying how the chiral phase transition proceeds through 
disoriented chiral condensates. Besides the assumption that the transition proceed slowly, no further assumption is 
introduced. In Section VII, we compute the effective action to the lowest one-loop order and derive the system of 
coupled equations that governs the evolution of the condensate fields, the equations which consist of the generalized 
relativistic Boltzmann equation, the self-consistent gap equation, and the equation of motion for the condensate fields. 
At the region where the curvature of the "potential" is negative, the effective action T develops imaginary part. It is 
shown that Im T can be dealt with by introducing the random- force fields that act on the system. The random- force 
fields cause negative correlation between the condensate fields. In Section VIII, we set up simple situations, for which 
approximate analytic calculations may be carried out. Section IX is devoted to concluding remarks. 



Let G(x, y) be (a component of) a generic two-point function. We introduce two kinds of spacetime coordinates: 



As mentioned ab ove , A-dependence of G{x, y) is assumed to be weak. Fourier transform of G(x, y) with respect to 



Throughout this paper, we use upper-case letters, P, etc., to denote a four-momentum vectors, P = (P M ) = (po,p), 
etc. The magnitude of three vector p will be written as p — |p|. As to the A-dependence of G(X;P), throughout 
this paper, we keep only up to a first-order-derivative term (gradient approximation), 



Remarks on notations and the gradient approximation 



x — y : relative or microscopic spacetime coordinates, 
x ~\~ y 

X = : center of mass or macroscopic spacetime coordinates. 




x — y reads (cf. (|l.f|) ) 
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G(X; P) ~ G(Y; P) + (X — Y) ■ d Y G(Y; P). (1.2) 

Functions of the type, i?(x, y; zq) also enter in Sections IV and VI. Here X (= (x + y)/2) and zq are macroscopic 
spacetime coordinates, i.e., i?(x, y; zo) depends on X and on zq only weakly. Fourier transform of -B(x, y; zq) on x — y 
reads 

B(x,y;z ) = / 7S5 e^-^X, z„; p). 

Whenever necessary, we expand B within the gradient approximation: 

B(X, z ; p) ~ B(Y, x ;p) + (z - x ) p) + (X - Y) • V Y B(Y, x ;p). (1.3) 

For the system of our concern, \dB/dxo\ < \B\/Lq, \dB/dx l \ < \B\jL % with L M as in Section I. 

In Sections VI - VIII, there comes in the condensate fields (p. Although the argument of <p is always macroscopic 
spacetime coordinates, we write <p(x) or <p(X), depending on the contexts. 



II. CLOSED-TIME-PATH FORMALISM 



Let us start with a generic Lagrangian density 

that governs the quantum-field systems. Here the superscripts label collectively the field type and internal degrees of 
freedom. Following standard procedure, we introduce an oriented closed-time path C (= C\ © C%) in a complex time 
plane, which goes from — oo to +00 (Ci) and then returns back from +00 to — 00 (C2). The time argument xq of the 
fields, (f>( a > (x)'s, is on the time path C. Along this closed time path, the classical action is defined as 

A = J dx J d 3 xC^ a) (x),(t> w (x),..^ 

= f £^(x),^\x),..) . (2.1) 

Define a path-ordered product of field operators along C: 

Tc^\ Xl ),...,^(x n )). 

Here Tc means to rearrange the field operators as follows: When Xjo of (j)^ aj \xj) is nearer to the end point of C, 
— 00 S C2, than Xio of (f>^ ai \xi), (j)^ aj \xj) is placed on the left of (f>^ ai \xi) with the sign due to the statistics of the 
relevant fields taken into account. Path-ordered Green function is defined as an average over ensemble prepared at 
an initial time t = — 00, which is characterized by a density matrix p: 

G = (-i)"- x Tr [lb (4> (ai \xi)-4 (an) {x n )) p 

ee {-iT-^Tc (4>^\x 1 )..4 (a - ) {x n ))). (2.2) 

In particular, the two-point function is defined as 

G^(x,y) = -i(TcU {a \x)4>« 3 \y))). 



Let us introduce || matrix notation G^ a ^ (x, y), whose component G I -" /3 ' ) (x, y) stands for G^ a ^ (x, y) with xq € C* 
and yo £ Cj\ 



G^\x,y) = -i(T (>)(x)<^(y) 

G[f\x,y) = - i a(^Hy)cf > ^(x)), 
G^\x,y)=-i{^\x)^\y)), 

G { 2 a 2 P \x,y) = -i(T(^)(x)0W(y))), (2.3) 
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where T (T) is a time- (an anti-time-) ordering symbol and a = — 1 when both <f>^ and <f>^ are fermionic field and 
a = +1 otherwise. In writing down (Pj), the "limit" < - Q - ) (xo € C\, x) = 0^ Q ^(a;o S C2, x) has been taken. Throughout 
this paper, we freely take this limit rl 

Let us break up the time path C into two parts: 



So 



dxo 









/•OO 




/ dx -f 


/ cfeo = / dxo 




/ da; 






./ C2 J ~oo 


10 £ ft 


' —OO 





(2.4) 



Then, in conformity with the matrix notation, we may rewrite the classical action in (2.1) as 

f 3 3 

cfcco d x C, 



A 



C = C[cj>[ a) {x)Af ] {x),. 



(2.5) 



where ^.j (2;) stands for the field whose time argument is on Cj (i = 1,2). The field j/)^ is called a type-1 or 
physical field while 4>2 IS called a type-2 field. We call t in ( pH| ) the hat-Lagrangian density |§,|l5|. Let us write 
(2.3) collectively 



G^(x,y)=-i{f (>)(*) *^(y))), 



(2.6) 



where 



ft a \x) 



4 a \x) 



t ^\y) = (4 m (y),4 m (y)). 



Four G l -" /3 ' ) (i,j = 1,2) in (|2.3|) are not mutually independent. As a matter of fact, it can readily be found from 



(2.3) that 



£(- 



Retarded, advanced, and correlation functions are related, in respective order, to G^'^ through 

G { f\x, y) = -i6(x ~ yo)( U [a) {x), ^(y)] ) 



G^'(x,y)-G[fHx,y), 
G^ \x,y)=i0( yo -Xo)(U^(x) )( j } ^(y)\ ) 



G[f\x,y)-G^>(x,y), 
-t(U a \x),cf>^(y)\ ) 

G^\x,y)+G^\x,y), 



(a/3) / 



G^\x,y) 
where [A, B\ a = AB + oBA (cr = ±). 



(2.7) 



III. COMPLEX SCALAR FIELD 

We take up self-interacting relativistic complex-scalar field theory with a conserved charge, whose Lagrangian 
density is 



Tn this respect, cf. the argument given below in conjunction with (7.49) 
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£ — £o + Ant ! 

Ant = -^(4^) 2 , (3.1) 

where the suffices "B" stand for bare quantity. Aiming at constructing perturbation theories, we go to an interaction 
picture. One can choose any value for the time t — U, when the interaction picture and the Heisenberg picture coincide. 
As in, e.g., pPJlq], we choose ti = — oo. For dealing with ultraviolet (UV) divergences, throughout this paper, we 
adopt the modified minimal-subtraction (MS) scheme [0, which is a variant of mass-independent renormalization 
schemes. We introduce UV-renormalized field <f>, mass m, and coupling constant A through 

cj) B = VZtfr, m 2 B = Z m m 2 , X B = Z\\. (3.2) 
Note that, although the same letter <f> is used, 4> here is the UV- renormalized field in the interaction picture, while 



>'s in Section II are the field in the Heisenberg picture. C in (3.1) may now be written as 



C = CM ,4>)+ Lint U\<f>)+L 



C rc {<t>\ <t>) = {Z- l)d^d»<t> - {Z m Z - \)m 2 ^^ - {Z X Z 2 - 1)~(^) 2 - (3.3) 

We shall treat Co as the nonperturbative part and Ant + C rc as the perturbative part. The UV renormalizability of 
the theory is demonstrated at the end of this section. 



The hat-Lagrangian density (cf. (2.5)) reads 

C = £(4,4> 1 )-£(4,<f> 2 ). 



From this with (3J3) we can read off the vertex matrices: Ant yields the (4>^(f>) 2 vertex matrix, 

iV 4 = -iXh (3.4) 

with T3 the third Pauli matrix and Ac yields 

iV 2 = -i [(Z - l)d 2 + (Z m Z - l)m 2 ] f 3 {4<j> vertex), 

iVl = -i(Z x Z 2 - l)Af 3 {{^4>) 2 vertex). (3.5) 

A. 2n-point Green function 

We assume that the density matrix p commutes with the charge operator, 

[p,Q] =o, 

Q = i ( d 3 x{tf<j> -$</>), (3.6) 



where <f> = d<p(x)/dxo, etc. (cf. below for undoing this assumption). Then, with obvious notation, nonvanishing 



path-ordered Green functions (cf. (2.2)) are 

G(l,...,n;^...,nO = i (-r(T c |0(l)...0H0t(l')...0t (n ' )e ^ c ^'| )5 

where A = Ant + C rc . Here, it is convenient to introduce a generating functional: 

Zc[J] = (T c exp (ij [C'(^(x), ( f>(x))+J*(x) ( p(x)+ ( pHx)J(x)] S j) (3.7) 



G 



where J(x) and J*(x) are c-number external-source functions. G is obtained from (3.7) through 

5 n Z c [J] 



G(l,...,n;l',...,n')=i 



5J*(l)...6J*(n)6J(V)...6J(n') 



j* = j = o 



(3.8) 



Hereafter we refer Zc [J] to as the CTP generating functional. It should be understood || that at the end of calculation 
we set J\{x) — J2(x), where Ji(x) = J(x) with i € Q (i = 1,2). 

As a generalization of the Wick theorem in vacuum theory, we obtain §,||Jl| 



T c exp(i [J*{x)<t>{x)+<tt{x)J{x)] 



c 



= p~ l fa fa J '( x ) A oc( x -y)J(y) 



= < : T c exp [i j [j*{x)<j>{x) +<f>\x)J{x)] ) :, 



(3.9) 



where the symbol : ... : indicates to take the normal ordering with respect to the creation and annihilation operators 
in vacuum theory. In (3^9), Aqc{x — y) takes the form, in matrix notation, 



A (x - y) 
Ao(P) 



d 4 P 



where 



(27r)4 e--^-)A (P) 
A F (P) 

e( Po )[A F (p) - a* f (p)] 



A F (P) 



p2 _ m 2 _|_ 



(-p )[A F (P) - A* F (P)} 
-A* F (P) 



(V = 0+) 



(3.10) 



is the Feynman propagator in vacuum theory. Simple manipulation yields 

(■Tee 1 fc^ J " '( x )<t>( x )+'t> i '( X ) J ( X )} . \ _ e iW c 

with 

Wc = '£— r (:T c \ [r(x)4>(x)+4t{x)J(x)] 



C 



n— 1 v ' 



W2n(xi, ...,x n ;yi, ...,y n )J*(xx)...J*(x n )J(yx)...J(y n ). 



Here the subscript 'c' denotes connected part and 

W 2n (xi, ■-,x n ;y 1 , ...,y n ) = i{-) n ( : (t>{xi)..4(x n )$ '(z/i). {Vn) ■ }< 



(3.11) 
(3.12) 

(3.13) 



represents initial correlation. In obtaining (3.12) from (3.11), we have used (3.6), which leads to 

(: (j)(x 1 )...(j)(x m )4> [ (yi)...^ (y n ) :) c = if m ^ n. 



If the condition (3.6) does not hold, the following argument in this paper should be generalized accordingly. Each 
of the 2n time arguments in (3.13) is on either C\ or C 2 - Then ( 3.13| ) consists of 2 2n components, all of which are 
identical. In the interaction representation adopting here, (<9 2 + rri 2 )<f)(x) = 0, from which we obtain 



(d 2 + m 2 )W 2n (xi, ...,x n ;yi, ...,y n ) = (d 2 + m 2 )W 2 n(x 1 , x n ; yi, y n ) = 



(j = 1, ...,n). 



The bare two-point function or propagator, A(x,y), is obtained from (3.7) with C — and (3.8) - (3.12) 



A{x,y)=G(x,y) 



A, 



£' = 
1 ±1 
±1 1 



A (x-y)+W 2 (x;y)A + , 



(3.14) 

(3.15) 
(3.16) 
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Thus we have learned that, in perturbation theory, the nonperturbative part consist s of t he bare two-point function 
iA(x, y) and the multi-point initial-correlation functions i(— )™ _1 W2n (n > 2). From ( 3.14 ), we see that W2n does not 
contribute to the on-shell amplitudes. Throughout this paper, we shall assume that WWs (n > 2) can be treated as 
(nonamputated) "perturbative parts" and, whenever necessary, we include "2n-point propagator" 

-i(-) n W2n(xu...,x n ;yx,...,y n ) (n > 2) (3.17) 

into the Feynman rules. We shall discuss, in the following subsection C, to what extent can one discard V\>2n (n >2). 



B. Bare two-point function 



We analyze in detail the bare two-point function (3.15), which is a building block of the perturbation theory: 

A(x,y) = -i{f (j>(x)4(y))) 



d 4 P 



{d 2 x + m 2 )A(x, y) = {dl + m 2 )A(x, y) = -f 3 S\x - y). 



Here Aq is as in ( 3.1C ) and 



= 02), ft = (4,4), 

W(x;y)=m(x;y)=-i(:<t>(x)4(y) :>c 
(Sg + m 2 )W(x; y) = (d 2 + m 2 )W(x; y) = 0. 



(3.18) 
(3.19) 



(3.20) 
(3.21) 



Let us find an explicit form for W(x; y) in terms of the density matrix p. For the purpose of later use, we consider a 
system enclosed in a rectangular parallelepiped. It is natural to take the sizes V (j = 1,2,3), introduced in Section 
I, as the lengths of edges. As discussed above, when we deal with microscopic reactions, IJ is regarded as large. For 
simplicity of presentation, we assume that L 1 = L 2 = L 3 (= L). For a set of single-particle wave functions, we employ 
the periodic boundary condition, on the basis of which (f>(x) may be expanded as 



= E 



(3.22) 



where V — L 3 , E p = \J p 2 + m 2 is the energy of <j> with momentum p, 

0,±1,±2,... (i = 1,2,3), 



2tt 

p = — n, n" 



L 



and a p (b p ) is an annihilation (a creation) operator of a particle (an antiparticle) of momentum p: 

Op|0) = hp|0) - 0. 

Substituting ( 3.22 ) and its hermitian conjugate <ft to (|3.20|) , we obtain, 



(3.23) 



(«t ap ) e - l ( p — Q-v) + ( 6q6 t ) e i(P-*-<3- 



y) 



where p° — E p and q° = E q . Taking the limit V (= L 3 ) — > oo and Fourier transforming on x — y, we obtain, after 
some manipulations, 



(3.24) 
(3.25) 



W(X; P)= J d 4 (x - y) e ip < x - y) W(x; y) 
F(P, Q) = S(2P ■ Q) 5{P 2 -m 2 + Q 2 /4) 



E ( T Po + r */ 2 ) 0(TPo - ^o/2) F r (rp - q/2, rp + q/2), 



T = ± 
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where 



F±(p,q) = (4(p)a±(q)). 



with 



a +(p) = 1™ 



I2E P V 



V^oa y (2?r) 



3 "P' 



o_(p) 



lim 



2E p V 

Or 



^oo V (2tt) 3 p 



(3.26) 



(3.27) 



In the case where the density matrix p is diagonal in Fock space, F±(p, q) oc S(jp — q). Substitution of thi s into 
( IH ) yields S(q), with which we have 5(2P • Q) -> 5(g )/(2|p |), so that .P(P, Q) cx (5 4 (Q). Thus, we see from ( |3~2l ) 
that W(X\ P) is X independent or W(x; y) (= W(ir. — y)) is translation invariant. The density matrix dealt with in, 
e.g., [Of is of this type. The case of equilibrium system is also of this type: 



F±(p,q) = 2E p 5(p - q)N ± (E p ), 

where -ZV+(-) is a particle (an antiparticle) number density. Equation ( [3.281 ) leads to 

W eq u{X;P) = W eQ u(P) 

= -2ki5(P 2 - to 2 ) [6(p Q )N + (E p ) + 9(-p )N^(E p )} , 

from which we have for the bare two-point function, 

A e?u (P) = Ao(P) - 2™<5(P 2 - to 2 ) [6(p )N + (E p ) + 9(- Po )N-(E p )} 



(3.28) 



(3.29) 



(3.30) 



a well-known form in equilibrium thermal field theory 

Let us turn back to the present out-of-equilibrium case As discussed in Section I, W(X;P) does not change 

appreciably in a single spacetime cell. Then, we see from (3.24) that, for jQ^I > 1/L 11 , T(P : Q) — 0, which means 



P±(±p-q/2,±p + q/2) ~ for q > 1/L. 
Here let us make following observation: 



(3.31) 



1) P is the momentum that is conjugate to the relative spacetime coordinates x — y, where x and y are in a single 

spacetime cell. 

2) In computing an amplitude in perturbation theory, W(X; P) is to be multiplied by some function P(P) and is to 

be integrated out over P M . 

The point 1) means that W(X;P) is meaningful only for \p°\ > 1/L° and p > 1/L. Then, the computation of the 
amplitude is reliable only when 

poo /*oo /*oo 

/ d 4 PW(X;P)F(P) ~ / dp d 3 pW{X;P)T(P). 

JO Jl/L° Jl/L 



Let us start from analyzing the "hard-p region," p > 0(1/L). Substituting (3.31) into (3.25), we see that P M in (3.25) 
satisfies 

pi = p 2 + to 2 + 0(1/L 2 ). 
To take care of the 0(1/L 2 ) term, we expand <5(P 2 — to 2 + Q 2 /4) in (3.25) in powers of Q 2 , 

2 Q 2 dS{P 2 - m 2 ^ 



5(P 2 - m + Q 2 /4) = 5{P 2 - to 2 ) + 



dP 2 



Substituting this back into ( p.2q ), we obtain for W(X;P) in ( 3.24 ), with obvious notation 

d 3 q 



W(X;P) 



-2m 



2 bo I 



ip-qX /po+«q-X 



Wpo + r|p • q|/2p )P r (rp - q/2; rp + q/2)] . 



(3.32) 



T=± 
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In con formit y with the gradient approximation we ar e adopting (cf. (1.2)), we neglect the term with d\ and with 
in (|3.32|) . Using 9(±p ± q Q /2) ~ 8(±p ) in ( |3.25| ), we finally obtain 



W(X; P) ~ -27r l( 5(P 2 - m 2 ) [0(;po)A + (A; P P , p) + 6(-p )N_ (X; E p> -p)] 



(3.33) 



where p = p/p and 



A±(A;P p ,±p) = y |^e-^ Xo ^ + ^ x P ± (±p-q/2;±p + q/2). 



Equation ( 3.33 ) is of the standard form, which is used in the literature. 

Now we turn to t he a nalysis of the region p ~ 1/L. It can readily be seen that the condition for W(X; P) to take 
the standard form ( 3.33 ) is 



F±(±p - p/2; ±p + p/2) ~ for q > \6i\ < 0(1/ L). 



(3.34) 



This condition is much stronger than (3.31). However, in most practical cases, requiring the weaker condition (3.31) 
is sufficient. This can be seen as follows. Let V be a typical scale(s) of the system under consideration. In the case 
of thermal-equilibrium system, V is the temperature of the system. As will be seen later, there emerges an effective 
mass Mi n d(X). In the case of m >> \[XP ', M inc i(X) is not much different from m and, for m < VXP, M in ^(X) is of 
0(V\V): 



M ind (X) = Max(m,VxV). 



(3.35) 



Most amplitudes, when computed in perturbation theory to be constructed below, are insensitive to the region 
\P»\ « Max(m,V\V) of W(X;P), and then the precise form of W(X; P) in this region is irrelevant. Here, it should 
be noted that m and \f\V are the scales that are characteristic to microscopic correlations. Then, from the setup in 
Section I, we assume that Ma x(m ; V~\V) > 0(1/ L). Thus, for dealing with the amplitudes of the above-mentioned 
type, requiring the condition ( 3.3l| ) is sufficient. 

Throughout in the sequel, for the piece W(X;P) of the propagator (cf. ( [3.18 ) and ( [3.24 )), we use ( 3.33j ) for all P M 
regions — oo < P M < +oo. As pointed out above, in computing some quantity, if infrared divergence arises at p = 0, 
we should reanalyze the quantity using the fundamental form, (3.24) with (3.25). 

After all this, we find for the Fourier transform of A(x, y): 



A(A;P)~ A (P)+W(X;P)A + , 



(3.36) 



with Aq(P) as in ( 3.10 ) and W(X;P) as in ( |3.33| ) . We refer the representation fl3.36j ) to as the FF representation, 
where F (F) stands for A F = (A ) n (A* F = -(A ) 22 ) in A (P), Eq. (pc|) . 

Infrared region in massless complex-scalar field theory, IZir, is defined in momentum space as IZir — {P p ; |P M | < 
O(VXP)}. In quantum electrodynamics (QED) and QCD, Kir = {P^ L ;0(g 2 T > ) < |P^| < O(gV)} with g the gauge 
coupling constant. It is well known that, in such theories as QED and QCD, the infrared region is full of rich structures 
and contains interesting physics JlJ,|l^,go| . Some quantities in QED or QCD is sensitive JlJ] to a yet "deeper" infrared 
region, W IR = {P^; 0(g 3 V) < |P^| < 0(g 2 V)}. If the inverse size of a spacetime cell, TZ ceU = {P^; \P^\ < l/L^}, is 
much smaller than the region IZir or 7Z' IR , then the above-mentioned interesting physics may be treated within our 
perturbative schemes. In the opposite case, the interesting physics is hidd en ou tside of our schemes. In more precise, 
the region IZir or TZ' IR is outside of the region, where the representation (3.36) is meaningful. 



C. Multi-point initial correlations 



In this subsection, we discuss to what extent V\?2n (n> 2), Eq. (|3 . 13[) , may be discarded. First of all, we recall that 
all the 2 2n components of W271 are identical. Then, from ( 
the retarded and advanced Green functions. This means that W2 does not appear in the two-point correlation in the 
linear response theory. With the aid of the standard formulae (cf. Sections 2 and 5 of |Q), one can also show that 
W2n does not appear in multi-point correlations in the nonlinear as well as the linear response theory. 

We restrict our concern to W4. The argument below may be generalized to the case of general W2n (n > 3). 
Straightforward computation using ( |3.2S| ) yields 



2.7) with Q3.1q), we see that W2 (= W) does not appear in 
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W4,(xi,x 2 ;yi,y 2 ) = i ^2 ^2 



pi, qi P2. q2 



n 



•'qiq2;piP2 



where pi = E P17 etc., and 

a 



qiq2;piP2 



= (at 



qi°q 2 a Pi a P2> ~ < a qi a PlX a q 2 a P2> ~ ( a q 2 °Pl ) ("4 fl P2 



(3.37) 



(3.38) 



Here we have assumed for simplicity that (a qi a q2 %>i%> 2 ) = (^ qi & q2 &pi^p 2 )j etc- The first exponential function in 
(3.37) may be written as, with obvious notation, 



where X — (x\ + x 2 + x 3 + x±)/A is the center-of-mass or the macroscopic spacetime coordinates and x'j and y'j arc the 
microscopic spacetime coordinates. Sim ilar decomposition can be made for the second exponential factor in ( [3.37| ). 
As setted up in Section I, W4 in ( |3.37 ) depends weakly on X M , which means that, for \p\ + p\ — — q®\ ^ 1/L° 
and/or |pi + p 2 - qi - q 2 | > 1/L, 



£-qiq 2 ;pip 2 — 0- 

Thus, the momentum conservation holds approximately in a single spacetime cell. Note that 



-'qiq2;piP2 



0(1) for 



3=1 



< 1/L° and 



X](Pi-qi) 

3=1 



< 1/L. 



(3.39) 



(3.40) 



Now we take the limit V — ► 00 and go to the momentum space: 

G A = -i J dx 1 dx 2 dy 1 dy 2 e^^ (P3 ' Xj ~ Q3 ' Vj) W 4 (x 1 ,x 2 ;y 1 ,y 2 ) 
(2tt) 4 V 2 



4 y/E pi E p2 E qi E q 



C qiq2; P lP2 II ( S ( P 3° ~ E Pi ) S ( q J° ~ E H )) 



3=1 



+C- Pl -p 2 ;- qi -q 2 II O^' + %3'0 + 

3 = 1 



(3.41) 



A transition probability or a rate of a microscopic reaction is related to an on-shell amplitude. (In the case of 
equilibrium thermal field theories, this relatio n is s ettled in |pl|.) Let us now analyze the on-shell (p| = E p ) amplitudes. 
As mentioned above and can be seen from (3.41), G4 per se does not contribute to the on-shell amplitudes. Then, 
W4 or G4 appears as a part(s) of an on-shell amplitude. Let us see the structure of a connected amplitude S4, which 
includes a single G4. We may write, with obvious notation, 



Si = (2tt) 4 



2 d 4 Pj d 4 Qj 



n 



L (2n) 4 {2*)* 



G 4 H6 4 [^p i +^(Q j -P j ; 

3=1 



(3.42) 



Here J2i Pi + Ylj=i(Qj ~ Pj) ^ s ^he sum of the momenta that flow in H. As mentioned above Pi + P 2 ~ Q\ + Q 2 , we 
may replace the delta function in (3.42) with <5 4 (^ ; Pi): 



S 4 = 



n 



d 4 Pj d 4 Qj 

L (27T)* (27T) 4 



G4 H. 



(3.43) 
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Let G 4 be the Fourier-transformed 2n-point Green function computed in perturbative scheme constructed below, 
which is the counterpart of G4. Note that G 4 , being of O(A) with A the coupling constant, takes the form 

G 4 = {2k) a 5\P 1 +P2-Q1- Q 2 ) G' A . 
The contribution ( 3.42Q with G 4 for G4 reads 

^ = (27r) 4 5 4 (£P*)£ 4) 



ol . d i P 1 d i P 2 d*Qi ~, 



(2tt) 4 (2tt) 4 (2tt) 



(3.44) 



Here we recall that the volume of the box, V — L 3 , can be regarded as large. Let us see the dependence of ( 3. 43] ) and 
( |3.44 ) on (the large) V. G' A and then 54 are independent of V . Owing to the properties ( 3.3S ) and ( 3.4C ), integration 
over c\2 in ( |3.43j ) yields a factor of l/V. Then, at this stage, £4 is 0(V) larger than S' 4 . Let us consider possible 
properties of (a] 11 al l2 ap 1 ap 2 ) . 

PI) Strong factorizability. There exist scale parameters 62 (l^l > 1/L), so that, for |pi — P2I ^ l^l or |qi — cj2 1 ^ |<?2 1, 

( a qi fl q2 a Pi a P2) - («qi a Pl)K2 a P2> + ( fl qi °P2 ) ( a q 2 «Pl ) ) ( 3 - 45 ) 

which leads to C qi q 2 ;p lP2 ~ 0. 

PI') Weak factorizability. There exist scale parameters 62 (|<$2| > so that, for |pi — P2I > l^l and\qi~ CJ2 1 ^ l^l, 

(|1|) holds. 



P2) Strong "short-range correlation". There exist scale parameters S3 (\Ss\ > 1/L), so that, for |pi — qi| > | ^3 1 or 
|Pi - qa| £ N, 

( a qi a q 2 a Pi a P2> - - ( 3 - 46 ) 

P2') Weak "short-range correlation". There exist scale parameters ^3 ( | ^3 1 > 1/L), so that, for |pi — qi| > l^j and 
|pi — q2| ^ \Ss\, ( 3.46| ) holds. From ( 3.39| ) and ( |3 .4C ) , this condition implies |p2 — qi | > l^l and |p 2 — q2 1 ^ l^l- 



Incidentally, PI) and P2) are not compatible with each other. In fact, when pi ~ qi and |pi— q 2 1 > > M qx ffel, l^l), 
P2) yields (a q i a q 2 a Pi a P2> - °- 0n the other hand > Pl)_or PI') leads to (|3.45|) . Owing to ( |3.34|) with (|3.26|) and 



(3.27), the second term on the right-hand side (RHS) of (3.45) vanishes but the first term does not, which contradicts 
( a qi a q2 a Pi a P2> - above. 
Let us enumerate the cases, where S4 and then W 4 may possibly be ignored. 



Property PI): 



S4 

si 



o 



\S2\ 6 V 
XT 3 



(3.47) 



Property P2): 



l± = n MX. 
SL A7« 



(3.48) 



Properties PI') and P2): 



S4 _ \S 3 \ 3 Mm(\6 2 \ 3 , \6 3 \ 3 )V 

S' A { xv 3 



(3.49) 
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Properties PI) and P2'): 



S 4 \6 2 \ 3 Mm(\6 2 \ 3 , \S 3 \ 3 )V 



(3.50) 



Properties PI') and P2' 



S4 



O 



\S 2 \ 3 \6 3 \ 3 )V 
XP 3 



(3.51) 



Let us inspect ( j3!47|) . If \5 2 \ & < 0(XP 3 /V) or \5 2 \ < C^A 1 / 6 y/V/L), G 4 or W 4 may be discarded. An exceptional 
case is where H in ( [[!§ ) is as large as cx A/(|5 2 | 6 ^) at Pi ~ P 2 ~ Qi ~ Q 2 with P x 2 = P| = Q\ = Q\ = m 2 . If 
\S 2 \ = 0(l/L), (Eyl becomes 



/ 1 



\xvr 3 



(3.52) 



and then W4 may generally be discarded. Similar inspection may be made for (3.48) - ( |3.5l| ). 
Let us analyze the special case where the density matrix p is diagonal in Fock space: 

fak \P\ n v) = ^kk'P„(«), ("k |Pl n k' ) = S ^'P n (b), 



where 



»/(*>) _ 



0, 1, 2, .... In this case, Green functions are translation invariant and independent of the ccntcr-of-mass 



or macroscopic spacetime coordinates (cf. above after (3.27)). Equation ( |3.37 ) becomes 



{2E p Vf 



[{a^al)-2(ala p f] 



x e 



^ e -iP-(x 1 +x 2 -yi~V2) _|_, 



where 'c.c' stands for the complex conjugate and (3.41) becomes 

dxidx 2 dyidy 2 e l ^- / ^= 1 ^ Pj x ' 



G 4 



Vi) 



W 4 



(3.53) 



(2tt) 12 <5 4 



(g«-<w)^^[<*C 



2«a Pl ) 2 ] 



x6 i (P 2 -P 1 )S i (Q 1 -P 1 ), 
where, for simplicity, we have assumed N p = 7V_ p . This leads to 



04 
S4 



O 



( 1 



\xvr s 



(3.54) 



(3.55) 



which is of the same order of magnitude as ( |3.52[ ). Thus, G4 or W4 may generally be ignored. Incidentally, for W 2n , 
we obtain, in place of (|||), S 2n /S 2n = 0(1 /(WV 3 )™- 1 ). 

In the sequel of this paper, we ignore throughout the multi-point correlations, W 2n (n> 2). 



D. UV renormalizability 

Here, a comment is made on the UV renor mali zabilit y. W e have completed the construction of th e building blocks 
of Feynman rules. The vertices are given by (3.4) and (3.5) while the propagator is given by (3.36). As discussed at 
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the end of Section IIIA, it should be stressed that, if the "2n-point propagator" (n > 2) W^n cannot b e disc arded, we 
should include it into the rules. We assume that, as in equilibrium thermal field theories, W(X; P) in ( 3.36| ) and WW 
(n > 2) in ( $1% are "soft." Here, by "soft" we mean that, in the UV limit |P^| -> oo, W{X; P) (e A(X; P)) tends to 
zero more quickly than the vacuum part A(°) (P) and W 2n tends to zero more quickly than the (leading contribution 
to the) 2n-point function. Then, genuine UV divergences emerge only from the vacuum part A^°^(P) of A(X;P). 
This means that Green functions are UV finite if the vacuum theory is renormalizable and has been renormalizcd, 
which is in fact the case for vacuum theory of complex-scalar fields. For more elaborate argument, we refer the reader 

to g|§. 



IV. QUASIPARTICLE REPRESENTATION 
A. Introduction of quasiparticle 

It is convenient to rewrite A(x, y), Eq. ( [3. IS ) or ( 3.36| ), using the so-called "Retarded/ Advanced basis" 
A(x, y) = A RA (x -y)- iF(x, y)A+, 
A RA (x -y) = A diag (x -y) + (A R (x - y) - A A (x - y)) ( ^ jj 

d 4 P 



F{x,y)=iW(x;y) + 



(2tt) 3 



-*P<*-v)0(-p o )S(P 2 -m 2 ). 



Here 



with 



A d iag(x -y) = diag (A R (x - y), -A A (x - y)) , 
A R(A) {x-y) = J ^0L e - lP -^A R{A) (P) 



Afl(A)(P) = ^2 2 _i_ ' / ^~ 

P z — m z ± ie(po)r) 



(4.1) 
(4.2) 

(4.3) 



(4.4) 
(4.5) 

(4.6) 



the retarded (a dvan ced) Green function. We refer (4.1) to as the R/A representation. When compared to the FF 
representation ( 3.36| ), the R/A representation makes the subsequent expressions compact and convenient for handling. 
A cost to be paid is that the R/A repres entat ion ge nerat es less intuitive expressions. 
In the case of equilibrium system (cf. (HJ|) and (|3~30|)), we have 



A e9 „(x - y) = A RA (x - y) 



d A P 



(2tt) 



- p <'-vh(po) M-po) + N(E P )} [A R (P) - A A (P)} A+, 



where a chemical potential (being conjugate to the charge) has been assumed to vanish, so that N+(E P ) — N-{E P ) 
(= N(Ep)). We note that, in momentum space, A equ may be diagonalized as 



A equ (P) 



fequ(P0,E p ) \i (p\f 1 + fequ(P0, E p ) fequ(po, E p ) \ 

+ fe q u(po,E p )) /Xd ^ [r) { 1 1 j 



where 



1 

1 1 + fequ{P0, E p ) 

f eq u(po,E p ) = e( Po )[6(- Po ) + N(E p )l 



(4.7) 



In configuration space, ([4.7|) reads 



A equ {x-y)= f AcI 3 ^bM( x -u)A<;» 9 (u-v 

J T = ± 



x - yo ) B^\v -y), 



(4.8) 
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where [A^j ] stands for a positive [negative] frequency part of Adi ag and 



1 J U(x) si*) + f£i(x) r * vy 



^(x)+/i 9 T 2(x) /M(x) 
5(x) 5(x) 



(4.9) 



C ( T ) M = ( d3p p'P x f 



with 



/ e ^(x)=e(r) y ^ e *P-*[0(- T )+7V(£ p )]. 

It is to be noted that the form for the matrices B^ T ,L, which leads to (4.7) with diagonal Adiag, is not unique and 
the choice here is the s imp lest one PJl5|. 

The representation (4.8) may be interpreted in terms of quasiparticle picture. In fact, (4.8) naturally leads us to 
introduce quasiparticle fields, 4>ra and (j> R A, through 

0(a:)= / d 3 u^Bi T) (x-u)41(u,x ) 
${y) = [ d 3 ^S(v,j/ )4 T) (v-y), 

J T = ± 



(4.10) 



where ^ra [4>ra] denotes positive frequency part of 4>ra [$ra] an d $ra [yra\ denotes negative frequency part. The 
transformation ( 4.1C ) is called thermal Bogoliubov transformation. The statistical average of T4>ra4>r A (cf. 
assumes the form 



{T(j) RA {x) 4> RA {y)) = iA d aig(x - y) 



(4.11) 



Thus the fields 4>ra and 4>ra describe well-defined propagating modes in the system, which we call retarded/advanced 
quasiparticles. It is to be noted that (j> RA is not a hermitian conjugate of 4>ra, which is a characteristic feature of the 
theory of this type 0Jl5 ] . 

With ( 4.1C| ) and ( [4.11 ), we can readily show that A equ {x — y) turns out to (4.8). From (4.E), we see that 



d 3 u B { L T} (x - u) f 3 B { ; ] (u - y) = f 3 <5(x - y) , 



(4.12) 



which shows that, under the transformation ( 4.1C| ), the canonical commutation relation is preserved: 



U x )> $ (y) 



S{x - yo) = 



4> RA (x), 4>R A {y) 

ir 3 S 4 (x - y), 



S(x a - y ) 



(4.13) 



where <fi (y) = dft (y)/dyo- Furthermore, the nonperturbative or bilinear part of the hat-Lagrangian density is form 
invariant, 



Co = -4 (d 2 + m')<t>x + {d 2 + m 



4> ^{d +m )4> 
4>RA f 3{d 2 + m 2 )(j)RA. 



(4.14) 



Let us turn back to the out-of-equilibrium case. As a generalization of (4.10), we assume that the fields, </> and </>' , 
are related to the retarded/advanced quasiparticle fields, 4>ra and $ra: through the transformation, 



x o)0fli(u, x ) 



kx)= I d 3 u^Bi T) (x,u 

J T = ± 

t( y )= / d 3 «^S(v )2/ o)5i T) (v,y;yo), 

J r=± 



(4.15) 
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where 



B[ T) (x, u;x ) 



B%\v,y;y ) 



£(x-u) /W(x,u;x ) 
<J(x-u) <5(x-u) + /( r )(x,u;a;o) 

^(v-y) + / (r) (v,y;yo) / (T) (v,y;yo) 
S(v-y) 5(v-y) 



(4.16) 



We assume that B^j) m (x, u; xq) depends on (x + u)/2 and on xq only weakly. Corresponding to (4.12), we have 



'L(R) 



d A u (x, u; x ) t 3 B%' (u, y; x ) = f 3 <5(x - y), 



(4.17) 



so that, as in (4.13), the transformation (4.15) preserves the canonical commutation relation. 

Incidentally, the so-called nonequilibrium thermo field dynamics 0,[l5| is formulated by taking (fOlj) as a starting 
hypothesis. 



B. Propagator 



Using ( 4.15 ) and ( 4.11 ), we obtain for the propagator A, 

= -if d 3 U d 3 t ;4 r) (x,u;xo)(f^(u,xo)^i t (v;y ))^ r) (v,y;y ) 

= f d\d 3 vj^ B l\^^o) A^ ag (u-v;x - y ) B^\w,y-,y ), 
J r=± 



which leads to (4.1) with 



F(x,y)=^F^(x,y) 

T = ± 

= i [ d 3 uJ2 [A^ ) (x-u;x -j/o)/ (r) (u,y;y ) 

J T=± 



-/ (t) ( x , u;a;o) A^ T) (u-y;a;o - y ) 



iW(x;y) 



d*P 
(2Vp 



,-iP-(x-y) n(_ 



p )S(P 2 ~m 2 ). 



The last line is the definition of F(x,y), Eq. ( |4.3|). E quation ( |4.19| ) is a determining equation for f( 7 
Straightforward manipulation using (1.2) in (4.18) yields 



where 



A(X;P) = I d 4 (x-y)e iP < x -^A(x,y) 

~ A(°)(X; P) + A^(X; P) (X = (x + y)/2), 
AC°)(X; P) = AiM(P) + A(°)'(X; P), 



* M j " ^ A R (P) - A A (P) -A A (P) 
A (0) '(X; P) = -27rte(po) / (e(po)) (X; p) 5(P 2 - to 2 ) i+ 

A«(X;P)4^^^- ( A fl( P) + A.(P))i + 
= -iP-9x/ (e(po)) (X;p) [A 2 f (P) + A 2 4 (P)] i+. 



(4.18) 



(4.19) 



(4.20) 



(4.21) 

(4.22) 
(4.23) 
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This is not the end of the story. We sho uld im pose t he co ndition (3.19), which leads to a constraint on f^ T \X] P). 
The most convenient way is to start from ( 3.21 ) with ( 4.19| ), 



{d 2 x +m 2 )F {T \x,y) = (d 2 + m 2 )F^(x,y)=Q. 
However, for the purpose of later use, we analyze ( |3.19] ) directly. From ( 4.18 ), we obtain 



(4.24) 



(d 2 + m 2 ) 



A(x,y) = (dl+m 2 ) I d 3 U d 3 «^B[ r) (x,u;x ) 

J T = ± 



xA di > a q ( u -v;xo-yo)B£ > {v-,y,y ) 



-t 3 S i (x-y) + j d 3 ud 3 v ^ 



8xq dxo 

{(X7i - V 2 U ) s[ T) (x, u; x )}] h B^iu, v; x Q ) f 3 A«(v, y; x ,y ), 



(4.25) 



where A^ + ) [A^ ">] is a positive [negative] frequency part of A. Comparing this result with ( |3.1S| ), we see that the 
second term should vanish, which reads, in the momentum space, 



2i 



d 4 P 
(2^ 



3 -iP-(x- 



x - y) J2e(r Po ) P-d x f^(X;p) AW(X;P) = 0, 



r=± 



where X — (x + y)/2. Since po and p are mutually independent, we obtain 

dx f {T \X; p) = p Vx/ (T) (1; P) = (r = ±), 



(4.26) 



(4.27) 



as it should be, since, interaction is absent here. Thus f^(X;p) is X independent and has vanishing gradient along 
the "flow lines." It can easily be shown that ( |4.27| ) guarantees also (d 2 + m 2 )F^ = 0, Eq. ( f4.24| ). Thanks to the 
relation (4.27), A' 1 ' in (4.23) vanishes, and then 

A(X; P) ~ A RA {P) - 2me( Po ) f^\X; p) 5(P 2 - m 2 ) A+ 

~ A RA (P) + [W(P; X) - 2nie(-p ) S(P 2 - m 2 )] A+. (4.28) 



|po| prescription 

Comparison between ( [4.28; ) and ( 3.33 ) yields 



f ( +\X;p) = N + (X;E p ,p) 

f(-\X;p) = -l-N-(X;E p ,-p), (4.29) 



where N + (N-) is the particle- (antiparticle)-number density function. 

In the case of equilibrium system with vanishing chemical potential (cf. ( |3.30 )) 



oE p /T _ i ■ 



fequ{E p ) always appears in combination with 5(P 2 — m 2 ) = 5{pj ) — E 2 ). Then, at first sight, it seems that no difference 

arises between fiqu{\po\) and feql(E p ). This is, however, not the case, since, in general, A(P) is to be multiplied 
by the functions that are singular at |po| = E p . The correct choice has been known p4j for some time now to be 

feql(\po\) — the "\p \ prescription." 

Let us turn back to the present out-of-equilibrium case. On the basis of the above observation, we assume that the 
\po\ prescription should be adopted: 
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/«(X; p) = A+(A; E p , p) -> A+(A;p , p) = /<+> (X; P) 
/ ( - } (A;p) = -1 - A_(A;£ p ,-p) -» -1 - A_ (X; -p , -p) = f^(X;P), 
/(A; P) = 9(p ) .f (+) (X; P) + 0(-p o ) / (_) (*; P). 



(4.30) 



In what follows, we simply refer / to as the number density (function) . An argument for the necessity of adopting the 
|po| prescription is given in Appendix A. 

It is to be noted that the translation into the |po| prescription is formally achieved by replacing (4.18) and ( 4. lGj ) 
with, in respective order, 



A(x,y)= Id ud v B L (x, u) A diag (u - v) B R (v, y) 



B L A 



L ■ ^diag 1 Br 



and 



B L (x,u) 
B R (v,y) 



S 4 (x — u) f(x, u) 

S 4 (x — u) 5 (x — u) + f(x,u) 

S 4 (v - y) + f(v,y) f(v,y) 
8 A {v - y) 5 A {v - y) 



The function f(x,y) here is an inverse Fourier transform of the number-density function f(X;P) in (4.3C): 

/(A; P)= ( d\x - y) e ip <*-^f(x, y) (X = (x + y)/2). 



Now (1.26) turns out to 



~ 2i J Wf e- iP <*-^ [P ■ d x f(X; P)\ A(X: P) = 0, 



and then, in place of ( 4.27 ), we have 



P-dxf(X;P) = 0. 



(4.31) 



(4.32) 



(4.33) 



This is a continuity equation for / along the "flow line" in a four-dimensional space. Equation (4.33) may be solved 



f(X-P) = F(X - (A - T m )p/ P0 ;P), 



(4.34) 



where T is an arbitrary function and Ti n is an (arbitrary) " initia l" time. Given the initial data f(Xo — Ti„,X;P) 
that characterizes the ensemble of the systems at Ao = Tj n , (4.34) fixes the form of T and we obtain 

/(A; P) = /(Ao = T^, X - ( X - T in )p/p ; P). 



The propagator A(A; P) with this /(A; P) still takes the form fl4.28|) : 

A(A; P) = A RA (P) + /(A; P) (A R (P) - A A (P)) A+. 



(4.35) 



In order to see the physical meaning of / to be the "number density," let us compute a statistical average of the 
current density: 



(*)) = 2 
_ 1 
~ 2 



(^(x) 0(a)) - {(j>{x) ${x)) 
(K - i^i2(y,x) + A 21 (y,x)] y=x 
i J [A 12 (x;P) + A 21 (x;P)], 



(4.36) 
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where d= d— d- Straightforward manipulation using dQ5| ) with ( p0| ) yields 

d 3 p 



<i°(z)> 



(2tt) 3 
d 3 p p 



where 



(2^) 3 £ p 



N + (x-p) = N+{x;P) 



N + (x;p)-N_(x;-p) 
N + (x;p)+N-(x;-p) 



(4.37) 



Po = E p 



N_(x;-p) = AT_0;-P) 



(4.38) 



Thus iV+(x; p) [iV_(x; — p)] is a number density of particles with momentum p [antiparticles with momentum — p]. It 
is to be noted tha t (|4.38p is valid in the gradient approximation, i.e., the terms with first derivative dN±(x; ±P)/dx 
do not appear in (|4.38[) . Also to be noted is that the argument x here is a macroscopic spacetime coordinates. 



C. Two perturbative schemes 



In this subsection, we introduce two perturbative schemes of calculating various quantities. The first scheme, which 
we c all the "b are-/ scheme," is the n aive one introduced above. In this scheme, perturbative calculati on go es with 
(3.4) and (3.5) for vertices and (4.35) for propagators. We designate the number-density function in (4.35) in this 
scheme as fs- 

The second scheme, which we call the "physical-/ scheme," is defined as follows. We redefine the fields, <f> and <jr , 
by ( 4.15 ) with ( 4.16 ), where is different from fg\ Then, the fields in this scheme are different from the fields 
in the bare-/ scheme. Changing the defini tion of fields leads to em ergen ce of a counter term in the Lagrangian. In 
order to extract it, we first note that ( 4.33 ) does not hold now and ( 4.25 ) may be written as 



/ d 3 ud 3 v^2D r (x;u : v)A^(v,y 1 x ,y ) =s -r 3 S 4 (x - y), 

J T = ± 



D T (x; u, v) = (fig + m 2 ) <5(x - u) S(u - v) - 

-{(V*-V*)Bf(x,u;z )} 
This means that, as a free hat-Lagrangian density, we should take 

C = 



dB^\x., u; xq) d 
dxo dxo 

T%Bfl ) (u ; v ; a; o)T3. 



d 3 u d 3 v M f (x; x )r 3 D T (x- u, v) 0« (v; x ) , 

T=± 



(4.39) 



(4.40) 



where cj>( r ^ — ((f)^^ , <j>2 ) anc ^ = > 4>2 ) with fy^ (<j>\ ■*) a positive (negative) frequency part of <pi, and 

a Hermitian conjugate of <^ . Inst ead of (4.25), starting with (9 2 + m 2 )A(x, y), we are led to t he sa me Cq as 
( 4.40 ) above, as it should be. Equation ( 4.40 ) tells us that, for compensating the difference between ( 4.4C ) and the 
original free hat-Lagrangian density, the counter term should be introduced in the hat-Lagrangian density: 



d 3 ud 3 v ^ (r)t (x;x o )f 3 



dBY>(x,u;x ) d 



dx 



dx 



-(V 2 - V 2 ) B[ T \ X , u; x )\ t 3 &%\u, v; x ) f 3 ^ T )(v; x ), 
which leads to a two-point vertex function, 



(4.41) 
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d 3 ? 



T=± 



dB^\x,u;x ) d . )(r , 



9x 



<9x 



(Vi-Vi)B^(x,u 5a :o) 



xf 3 B^ T) (u,y;xo)f3(5(a:o - y ) 



T = ± 
d 4 P 



(2ny 



»5^£-(*-v",)/<'W,..> 

■ e- <p ^-^ ^ 0(r Po ) [P • d x /M(X; p)l i_. 

r=± 



<5(x - y ) 



Here X = (x + y)/2 and A_ is denned as in ( 3.16 ). Going to the |po| prescription, we have (cf. ( |4.3C| )), 

d 4 P 



iVc{x,y) 



(27T) 



_ e -*P-(*-y) [p.d x f(X;P)]A- 



(4.42) 



The form for f(X; P) will be determined later so as to be the physical "number density." If we impose (4.33), (4.42) 
vanishes as it should be. 

The "free form" of the propagator ( 4.28) ) should be replaced by (4.2C) - (4.23) with the \po\ prescription: 



A(x,y) 



d 4 P 



■ e 



-iP-(x-y) 



A(X;P), 



(2tt)* 

A(X; P) = A ( °) (X;P) + A« (X; P) 
A<°> (X; P) = B L (X; P) A diag (P) B R (X; P) = A RA {P) + A (0) ' (X- P), 



B L (X;P) 



1 f(X;P) 
1 1 + /(X;P) 



B R (X;P) 



l + .f(X;P) f(X;P) 
1 1 



where A RA (P) is as in ( f4.2l| ) and (cf. ( |4.30D ) 

AW (X; P) = /(X; P) (A R (P) - A A (P)) A+, 

= -iP ■ d x f(X; P) [A|(P) + Ai(P)] A + . 



(4.43) 

(4.44) 
(4.45) 

(4.46) 



(4.47) 
(4.48) 
(4.49) 



Let us discuss two perturbative schemes in analogy with mass renormalization in vacuum theory. We take the 
Lagrangian density 



C — Co + Ci 



,1 a\, 



Co = d»<f> - m 2 B , 

The "bare scheme" in vacuum theory takes Co as a nonperturbative part. Then, the propagator reads, 

1 



A(P) = 



p2 



m B + i- 7 ! 



whic h corr esp onds to ( 4.2S ) in the bare-/ scheme here. The "renormalized scheme" in vacuum theory adopts (cf. 
(O and Q) 



C = d^(j} - m' 



(4.50) 



1/2 

as the nonperturbative part. Here m = Z m tub is a renormalized mass. Equation (4.50) leads to 

1 



A(P) = 



P 2 — to 2 + irj ' 



for the propagator, which corresponds to (4.43) - ( 4.49| ) in the physical-/ scheme here. To compensate the difference 
Cq — C'q, a counter term should be introduced into the Lagrangian density, C c = — (Z m Z— 1) m 2 ^4>, which corresponds 
to ([4.41|) here. 
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To summarize, the bare-/ scheme is constructed in terms of original bare number density Jb(X;P) [respect, bare 
mass tub]- On the other hand, the physical-/ scheme is constructed in terms of physical number density f(X;P) 
[respect, renormalized mass to]. Both perturbative schemes are equivalent. The first scheme starts with the "bare 
quantity" and the "renormalization" is done at the end, while in the second scheme, the "renormalization" is done at 
the beginning by introducing the counter hat-Lagrangian L c [respect. C c \. It is worth recalling that the renormalized 
mass to is defined so that to is free from UV divergences. However, there is arbitrariness in the definition of the "finite 
part" of to. to is determined by imposing some condition. (In this paper, we are adopting the MS scheme.) Thus, 
to is determined order by order in perturbation series. As will be seen in the next section, this also applies to the 
present physical-/ scheme. Namely, there is arbitrariness in defining /. In Section VC, we shall impose the con dition 
for determining / so as to be the physical number density, under which P ■ dxf(X; P) in A^(X; P), Eq. ( I.49| ), and 
then also f(X; P) turn out to be determined order by order in perturbation series. 

V. INTERACTING FIELD 

Interactions among fields give rise to reactions taking place in the system, which, in turn, causes a nontrivial change 
in number density of quasiparticles. In this section, we analyze full propagator in two perturbative schemes presented 
at the last section and demonstrate their equivalence. 

A. Self-energy part 

The self-energy part £ takes the form 

£(x, y) = t int (x, y) + t c (x, y), 

t c {x,y) = -V c (x,y) = -V c (x,y)A_, (5.1) 



as m 



where £mt is a loop-diagram contribution, £ c comes from the counter hat-Lagrangian density £ c in ( 4.41 ), and V c is 

(ED- 



It can be shown that (cf. Appendix B) 



£ E« - ]T = 0, (5.2) 



with which we obtain, in the configuration space, 



^B R -^-B L = ( ), (5.3) 



where P.r(z) is as in (4.32) and 

£fl = Ell + £l2 = (Emt)ll + (£mt )l2, 

£.4 = — (E22 + E12) = — {(Ej n t)22 + (Ej„t)l2} (= Efl), 

E o// - Eia + £12 • / - E 2 i • / + E A • / - / • E A + V. (5.4) 

Efl and Ea are called the retarded and advanced self-energy parts, respectively. 
Going to momentum space, we obtain, after some manipulations, 

E o// (A; P) ~ -i {f{X; P),P 2 -m 2 - ReZ R {X; P) } + ifW (X; P), (5.5) 

where 

Li(X-P) n(X-PV. - 9A(X;P)dB(X;P) dA(X;P) dB(X-P) 

{MX,P),B(X,P)}=—^ — — ^— , (5.6) 

if W(X; P) = (1 + f(X- P))£i 2 (X; P) - f(X; P)E 21 (JC; P). (5.7) 



In deriving (5.5), use has been made of the relation Ea(A; P) = \T,r(X; P)]*, a proof of which is given in Appendix 
B. 
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B. Self-energy- part resummed propagator 



In this subsection, we compute a self-energy-part resummed full propagator trough the Dyson equation, 



G(x,y) = A(x,y) + 



A-E-G 



We w ould like to obtain G(A; P) (A = (x + y) /2) , the Fourier transform of G(x, y) on x — y. Using the transformation 
( 4.31 ), we obtain 



G(x,y) 



B L 1 ■ G ■ B R 



(x,y) 

= A diag (x,y) + A diag -t-G (x,y). 
Fourier transforming on x — y using ( |C2| ) in Appendix C, we obtain 
G(A; P) ~ A diag (P) + A diag (P)t(X; P)G(A; P) 



+ - 



dA diag <9E - ~ ffi 3G dA diag Z dG 

Lr — l\ d i ag - 



dP^ dX„ 



' dX^ dP^ 



dP» dX„ 



(5.8) 



where (cf. (|5.3|) ) 



S — ^diag + E D // ( q q 

Sdia 9 = diag (E fl , -E^) • 



(5.9) 



The terms in the square brackets in ( |5.8| ) involve ^^-derivatives and are small when compared to the first term on 
the RHS. 

Let us solve (5.8) iteratively. In doing so, we assume that the leading terms on the RHS of are A d i ag + 

Adiag^diagG_ and E's in the square brackets may be approximated as E ~ Yi d i ag . This assumption will be justified a 
posteriori (cf. the next subsection). Straightforward manipulation yields 



where 



G(A;P) =G diag {X;P)-Z off (X;P)G R (X;P)G A (X;P) ( " Q 



G diag (X; P) = diag (G R (X; P), -G a {X; P)) , 

1 

Gria) (A; P) 



1 



(5.10) 



P*-m?-Y. R{A) {X;P)±ie(p )r 1 



(5.11) 



It is to be noted that ( p.lC| ) is exact to th e de rivative approximation, i.e., no term including the first derivative (with 
respect to A M ) arises. Then, using again (CS) in Appendix C, we obtain 



G(A;P) = G (0) (A;P) + G (1) (A;P), 
G< >(A; P) = B L (X; P)G dlag (X; P)B R (X; P) 
= G RA (X;P) + G^'(X;P), 



where 



Gra(X-P) 



G R {X-P) 
G R (X; P) - G A (X; P) -G a (X; P) I ' 



G (0) ' (A; P) = /(A; P) (G R (X; P) - G A (X; P)) A, 
G«(A;P) = -Z off (X;P)G R (X;P)G A (X;P)A^ 



(5.12) 
(5.13) 

(5.14) 

(5.15) 
(5.16) 



-- [{/(A; P), P 2 - m 2 - E«(A; P)} G R (X; P) 
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{f{X;P), P 2 



E A {X;P)}G 2 A {X-P)} i H 



^E o// (X; P) (G R (X; P) - G A (X; P)f A + 
1 

~2 

-if P) P) + Gi(X; P)) i + . 



{/(X; P), ImZ R (X; P)} (G 2 R (X; P) - G 2 A (X; P)) i+ 



(5.17) 



In narrow-width approximation, ImY, R = —ImT, A — > — e(po)0 + , the first term on the RHS of ( 5.17 ) develops pinch 
singularity, while the rest turns out to th e wel l-defined distributions. 

For clarifying the physical meaning of ( 5.17 ), we compute the contribution to the ph ysical number density through 
analyzing the contrib ution to the statistical average of the current dens ity (i ^ 1 ), Eq. (Of]). The second and third 
terms on the RHS of ( 5.17 ) leads to perturbative corrections to (j M ) in ( 4.37 ) due to quantum and medium effects, 
while the first term yields a large contribution, 



J^P" Z off (x; P) (G R (x; P) - G A (x; P)f , 



(5.18) 



which diverges in the narrow-width approximation. (Note again that the argument x here is a macroscopic spacetime 
coordinates.) In the next subsection, we inspect this large contribution more closely in the bare-/ and physical-/ 
schemes. 



C. Two perturbative schemes revisited 



1. Bare-f scheme 



Let us estimate ( 5. 18] ). All the formulae displayed above may be used with / — > /s and with imposition of ( 4.33| ). 
Important region of integration is where Re(G] i 1 (x; P)) = P 2 — m 2 — ReT. R (x; P) ~ 0. We define "on the mass-shell" 
Po = ±u±(X;±p) (= ±lu±) through 



[P 2 -m 2 -ReE R (X;P))] 



0. 



(5.19) 



We make following approximations; 

p"e 0// ~ 0( PO ) k, P )" si 6 ;; e)(+) + o(-po) (- W -, P )" ^; r f e){ -\ 



[G R {x;P)-G A {x-P)] 



0(t Po ) ImSg 



% [2u T Zr X (po - tlu t )Y + (JmS^) 2 ' 



where 



,(6are)(±) _ 
J off 



E o// (X;po = ±W±,p), ^ R = Y, R (X;p a = ±uj±,p), 



Z+ 1 = 1 T 



<9Pe£, 



po = ±u>± 



(5.20) 

(5.21) 
(5.22) 



It is convenient to carry out the wave- function renormalization such that Z± — 1, which we assume to be done for a 
while. Using (5.2C) in (5.18), we obtain 



1 f d 3 p 



2 J (2tt) 3 



(l,p/^+) A 



-iE 



(bare)(+) 
off 



\ImE 



(+)i 



(1,-p/^-) a 



(bare)( — ) 
off 



I/toE 



(-)i 



(5.23) 



In deriving this, we h ave a ssumed that |im£^| << ui±, which is the case in equilibrium case. As in (JD§) - ( pi ), 
we can read off from ( [5.23; ) the contribution to the physical number density: 
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5N±(x,±p) 



, v (6are)(±) 
1 ~ llj off 



(±) 



= 2 T± 



,-t(p) 



-ir 



2w± 



{/, tfeEfl} 



(5.24) 



Here r+ = 2cj + /|/mE^" ) | [r_ = 2a;_/|ImE '|] is a characteristic time during which a particle with momentum p 
[an antiparticle with momentum — p] damps: 



r A 3 

G fl (X ;2 ; -y ;p)= J — e^*"^ G R (x, y) 



(2tt) 3 

-Oo-J/o)/t+ 
-(xo-Vo)/t- 



for a particle mode, 

for an antiparticle mode. 



In (5.24), 



-.0) _ 



2ui± 



[(1 + iV±)E 12(21) - iV±£ 



J 21(12). 



Po = ±0>± 



(5.25) 



which comes from T^ p ' € —iYi tj in (5.7), is a net production rate of a particle (an antiparticle) of momentum p 



(-p). In fact, Y { l\ Eq. (^25|) , is a difference between the production rate and the decay rate, so that T±^ is the net 



production rate. In the case of equilib rium system, T^± — (detailed balance formula). 

To disclose the physical meaning of (5.24), for the tim e bein g, we put aside the s econ d term in the square brackets, 
whose physical interpretation will become clear later at ( 5.32| ). Then, the RHS of (5.24) is the change in the physical 
number density, during the time interval r±/2, due the the net production rate. In Section I, we have introduced 
spacetime cells, whose size is L M (/i = 0, 1, 2,3). It is quite natural t o tak e t± — Lq, where Lq is the size of a time 



direction of the spacetime cell (including the spacetime point in ( 5.24j) ). [Strictly speaking, i n gen eral, t + ^ t. 



However, in the present crude argument, we ignore this difference.] It is interesting to note that (5.24) is half of the 
net production probability, due to the reaction, during the time interval t± = Lq. [In this respect, cf. [K5J .] 

Noticing that we ar e concerned about the particle mode with momentum p [antiparticle mode with momentum 
— p], we see that ( 5.24 ) leads to the relations for the physical number density, 



Sibare) / 
Iff >{x- 1 u + ,p) 

iY^ff e \x-,-u-,y) 
2uj- 



N+{x + L , x + L p/u + ; p) - N+(x , x; p) ~ L Q - 



N^(x + L ,x - L p/^-; -p) - AT_(x ,x; -p) ~ L Q - 



from which we obtain 



dN ± (x;±p) 
dx Q 



± 



p dN± 

LU± (9x 



-iE 



(bare) 
off 



(x;w±,p) 



2ui± 



(5.26) 



Physical interpretation of these relations will be discussed below. 

In calculating so me qu antity in the bare-/ scheme, one carries out perturbative calculation using the propagator 
A(X;P) given in ( 4.35| ). The calculated quantity is written in terms of the bare number density /b(X;P) or 
Nb±(X; E p , ±p). Using the solution to ( 5.26 ), one can rewrite the quantity in terms of the "physical number density" 
N±(X;±p). 



2. Physical-f scheme 



As defined in Section IVC, we change the definition of the number-density function /. Then, in this physical-/' 
scheme, the number density function / does not obey (4.33). It is convenient to define /, such that the number 
densities 
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N + (X;E p ,p) 
N_(X;E p ,-p) 



f(X;p Q = E p ,p), 
-l-f(X; Po = -E p ,p), 



are as close as the physical number densities N±(X; ±p). As seen above, the first term on the RHS of ( 5.17 ), which 
is proportional to T, ff(X;P), yields a large contribution to the physical number density. Thus, we demand that 
f(X;P) changes in time Xq so that E Q //(X; P) vanishes: 



^ off (X;P)^{f(X;P),P 2 



ReZ R (X;P)}-tW(X;P) 



0. 



(5.27) 



Given an initial data at the "initial" time Xq — Ti„, f(Xo = Ti„,X;P), ( 5.27 ) serves as determining f(Xo,X;P) 
at later time Xq. Note that i?eE# and r( p ) are evaluated order by order in perturbation series. Then, through 
( |5.27 ), f(X;P) is also determined order by or der. (In this relation, see the argument at Section IVC for the mass 
renormalization in vacuum theory.) We refer (5.27) to as the "renormalization condition" for the number density, 
whic h re places the condition (4.33) in the bare scheme. In the present scheme, transformed self-energy part E, 
Eq. (5J5), is diagonal: 



BrZB l 



Z r 

-E A 



(5.28) 



It is to be noted that, with the condition (5.27), G(X\P) in ( 5.12 ) - ( 5.17 ) is free from singular contribution in 
the narrow- width limit (cf. above after (5.17)) and the difference between G and the bare propagator A leads to a 
perturbative correction to the physical number density. 

In the sequel of this paper, we adopt the physical-/ scheme. 

Let us inspect the physical implication of the "renormalization condition" ( [5.271 ) : 



2P ■ d x f(X; P) - {f(X; P), Re^ R (X- P)} = f <*> (X; P). 



(5.29) 



Here let us recall the definition (5.19) of the "mass-shell" and find the form of (5.29) on the mass-shell. From ( |5.1E| ), 
we have 



2p - 

8ReT, R 
OX 



dRe^j 
dpo 



v± - 2p - 



<9i?eE 



R 



dp 



dw ± ( dReZ R 



ax 



dpo 



po = ±ui±(X; ±p) 



±w±(X;±p) 



= 0, 



(5.30) 



where v± = ±duj±(X; ±p)/dp is the velocity of the ± mode. We also recall that / = 9(po)N+ — 9(— po)(l + N-), 
where 7V + = N + (X;po = w+,p) [N- = N-(X;— po = u-,—p)] is the number density of a (quasi) particle [an 
anti (quasi) particle] with momentum p [— p]. It is straightforward to show that, on the mass-shell, po = ±lu±(X; ±p) 
(= ±uj±), ( [5.29 ) becomes 



dN± 
~dX^ 



v± • V X N± ± 



8lo± dN ± 



dX^ dPf* 



Po = ±OI± 

dN±(X;tu±(X;±p),±p) duj ± dN± 



(p) 



Z±T ± 



dX 



du± dN± 

~dp ' ~dK ~ ~olt ' ~dp~ 



(5.31) 



where Z± is the wave-function renormalization factor (cf. (5.22)) and, as defined in ( 5.25| ), T± is the net production 
rate of a particle (an antiparticle) with momentum p [— p]. If we employ the wave-function renormalization condition 
Z± = 1, the RHS turns out to the net production rate r^. The qua siparticle distribution function N± — N±(X;po — 
uj± , p) here is nothing but the relativistic Wigner function and ( p\3l| ) is the generalized kinetic or Boltzmann equation 
(cf. ]2q]) for the relativistic complex-scalar-field system. 

Here we make a comment on a role of the second term on the left-hand side (LHS) of ( |5.29| ). The second term 
represents the effect due to the change of "mass;" 
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{f(X;P), ReZ R (X;P)} = 



df(X; P) dReZ R (X; P) 8f(X; P) dReZ R (X; P) 



dX u 



dPv 



dpn 



dX b 



(5.32) 



The first term on the RHS togeth er w ith the first term on the LHS of ( 5.29| ) yields the term being proportional to 
the first two terms on the LHS of (5.31), while the second term on the RHS yields the term being proportional to the 
third term on the LHS of J5 3l| ). 



One can show using ( |5.3C| ) that (5.26) is the same in form as (5.31). [Recall that we have s et Z± = 1 in (5.26).] 
This observation gives a support for the qualitative argument made above in conjunction with (|5.26j). 



D. Resummation of the lowest-order E 

To lowest order in A (cf. ( [3.3| )), a tadpole diagram contributes to £,- n t: 
±int(X) = M 2 (X)t 3 , 

M 2 (X)^i\ / ^A^>(X;P) + (Z m Z- 



(5.33) 



l)m 2 



p2 _ m 2 _|_ ^ 

-2m {6(po)N + {X;p , p) + 9(- Pa )N_(X; - Po , -p)}] S(P 2 
+ (Z m Z-l)m 2 . 



(5.34) 



Here use has been made of the dimensional regularization: D is the spacetime dimension and e = 4 — D works as 
a regulator for the UV divergence. For the purpose of later use, we have adopted here the FF representation (cf. 
( ^36| )) for A$, (fyij) with ( |4~45| ) and ( |4~47| ). N ± in is as in (^30|) . Note that A$ in ( |Q| ) with ( ^49| ) leads 

to a contribution, which is of higher order in A. It should be kept in mind that, in dealing with 0(A 2 ) contributions, 
(5.34) with for A^, as well as the two-loop contribution, should be taken into account. 



Standard manipulation yields 

M 2 (X) 



Ml 



X 



-^induced ("^0 



= A 



(4tt) 2 
d 3 p 



, to 
1 - In 



1 



(2tt) 3 2E P 



where 



D 



[N + (X;E p ,p) + N_(X;E p ,-f>)} 



(5.35) 
(5.36) 

(5.37) 



Here Mq is the squared self mass in vacuum theory, 7^ is the Euler constant, and is an arbitrary scale parameter. 
As an UV renormalization scheme in vacuum theory, we are adopting the MS scheme. Thus, at the present lowest 
order, we have 

A 2 



(4tt) 2 



and then 



M 2 = 



(4tt) ; 



-to 2 In 



TO 



(5.38) 



If we choose fid = m/^/e, we have Mq = 0. 

In passing, for a locally thermal equilibrium system, we have 

N+(X;E p ,p) = N-{X;E p ,-p) 



1 

3 E P /T(X) 
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where T(X) is the local temperature of the system. In this case, we obtain 

A 



M induced( X ) = ^2 m ~ 



. dU e mu/T(X) _ I ■ 



In the high-temperature regime, T(X) >> to, we have 



M induced( x ) — j^T 2 ( x ) _ n T(X) 



(T(X) » to), 



while, at low temperature, we have 



M^.,,X)^(Iffl) S,2 e-.»/^) 



3T(X)\ 



8 TO 



(T(X) « to). 



(5.39) 



(5.40) 



Let us construct a S-resummed propagator. To the present lowest order, £12 = £21 = and then, from (5.7), 
vanishes. Then, the renormalization condition (5.29), £ Q // = 0, reads 



whose approximate solution is given by 

f(X;P) ~ .P(X;P 2 - M 2 (X),p) ~ G(X;P). 
Here, T or G is an arbitrary function and 

dM 2 (X) {X -T m ) 



(5.41) 



(5.42) 



Po=Po 
Pi=Pi 



dX° 2p Q 
dM 2 (X) (X -T m ) 



dX i 



2po 



(no summation over i), 



Y _ Y Pi (Y rp v , Po9A/ 2 (X)/^- K c)Af 2 (X)/ax" 2 

— po 4(p r 



(5.43) 



Given an initial data f(X = T m ,X;P) = G(X;P), (|5.42|) gives /(X ,X;P) at later time X to 0(A). 

Since 7ro£ r = — 0, & 1 ' in ( 5.17 ) vanishes due to the condition ( 5.41 ). Then, the £-resummed propagator 
G(A~; P) reads 



G(X;P) ~ G (0) (a~;P) 

= G fiA + /(A-;P) (Gfl-G A )i+, 



(5.44) 



where 



G 



RA 



G 



R{A) — ^R(A) 



G R 
Gr — Ga —Ga 

Gria) [X] P) 



pi _ m 2 _ m 2 (X) ± ie(po)T?" 



(5.45) 



E. Changing the mass parameter 



The analysis in the last subsection shows that an efficient perturbative scheme is obtained as follows. Rewrite the 
Lagrangian density (3.3) by introducing a new "mass function" M(x), which is assumed to depend weakly on x: 
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where 



C = d^{x) d^(x) - M 2 (x) <pHx) 4>(x), 

C-int =~ ((p f (x) (/)(x)) 2 , 

C rc = X (x) ^(x) 4>(x) + (Z- l)d^(x) d^(x) - (Z m Z - l)(M 2 (x) - X (x)) <P\x) <f>(x) 
-\{Z x Z 2 -l)(4?(x)4>{x))\ 



x(x) = M (x) — to 



(5.46) 



(5.47) 



Note that M 2 here is not M 2 in (5.34). As discussed in pOj , perturbation theory based on (5.46) is UV-renormalizablc. 
We give an simple argument for this in Appendix D. 



It should be emphasized pQ| that the perturbation theory based on (5.46) is defined by taking Cg as a nonpertur- 
bative or "free" part and £i n t + £rc as perturbative part. Then, x{ x )4^ 4 1 m •Crc is a part of the perturbative terms. 
Thus, in traditional perturbation expansion, we assign 



M 2 {x) = 0(A°) , (= M 2 (x) - to 2 ) = 0(A), 

and, in loop or 6 expansion (cf. the next section), we assign 

M 2 {x) = 0{6°) , x (a:) (= M 2 {x) - to 2 ) = 0(6), 



(5.48) 



(5.49) 



It is o bvious from the argument in the above subsection that, in the bare-/ scheme in perturbation theory based 
on (|5.46|), the bare propagator is given by (5.44), 



where 



G(X; P) = G RA + f(X; P) (G R - G A ) A+, 



Gr(A) = G R (A){X; P) 



P 2 - M 2 (X)±ie(p )rj' 



(5.50) 



(5.51) 



The number-density function f(X;P) in (5.50) obeys (5.41). 

This result may al so be derived through the same procedure as in Section IV. This time we start with the |po|- 
prescription-formula (4.31) from the beginning, 



G(x,y) = / d 4 ud 4 v B L (x,u)G d iag(u,v) B R (v,y). 



Here Gdiag = diag(Gfl, — Ga) with G R ( A )(x, y) a solution to 

(a 2 + M 2 {x)) G R[A) {x,y) = (0 2 + M 2 (y)) G R(A) (x,y) = -S 4 (x - y), 



(5.52) 



(5.53) 



with retarded (advanced) boundary condition. It can e asily be shown that the solution to ( 5.53| ) is ( 5.51 ). Note that 
no first X^-derivative term appears there. In place of (4.25), we have 



(d 2 x + M 2 (x))G(x,y) = -T 3 6 4 (x-y) + J d 4 ud 4 v [{(d 2 - d 2 u ) 

+M 2 {x) - M 2 {u)} B L (x, u)} t 3 B R (u, v) f 3 G(v, y), 



(5.54) 



and in place of (4.41), we have 

to = - J d 4 ud 4 vft(x)f 3 [{(dl-dl) + M 2 (x)-M 2 (u)}B L {x,u)\ t 3 B r {u,v)t 3 4>(v). 
This yields a two-point vertex function 
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iV c (X;P) = 2P-d x f(X;P) 

= {f(X;P),P 2 -M 2 (X)} A 



dM 2 (X) df(X;P) 



dX^ dP» 



A. 



(5.55) 



Since we are to obtain the bare propagator, Sj„t = f ^ = 0. Then, the "renormalization condition" ( 5.27| ) turns out 
to 



{f(X;P), P 2 -M 2 (X)}=0, 



(5.56) 



with which we have V C (X;P) = 0. Equation ( |5.56 ) is nothing but ( |5.41 ). Using this condition, we can derive (5.5C) 
from (5.52), which is the bare propagator in the bare-/ scheme. The physical-/ scheme may be constructed in a 
similar manner as in Section IVC. 
In this theory, £j nt to O(X) reads 



^int(X) = f 3 



A . M 2 (X) , , N N , 

M 2 (X) In — ^ + M^ duoed (X) - M 2 (X) + m 2 



(47T) 



where M 2 nduced (X) is as in (|5.36[) with M 2 (X) for m 2 and — f 3 [M 2 (X) — m 2 ] has come from 4> € £ r c, Eq. ( |5.4t| ). 

Up to this point, the function M 2 (X) (> 0) is arbitrary. An efficient perturbation theory is obtained by demanding, 
e.g., t mt {X) = 0; 



M 2 {X) 



1 - 



A ^ M 2 (X) 



(4tt/ 



m 2 + M 2 



induced 



(X). 



(5.57) 



The scheme with this M(X) may be called "presummation scheme." Equation (5.57) is the gap equation, which 
serves as determining M 2 (X) in terms of fid, which is still an arbitrary param eter . 

It is worth mentioning that, in the case where the condition [p, Q] = 0, Eq. (3.6), is not fulfilled, we may generalize 
(5.46), although not necessary, by changing Cq to 



and changing C rc accordingly. Here the superscript "(+)" ["(—)] refers to the positive [negative] frequency part, and 
then M + (x) [Af_(x)] is an effective mass of a quasiparticle [anti-quasiparticle] . 



VI. THE O(N) LINEAR a MODEL 



A. Preliminary 



In this section, we deal with a fermionless O(N) linear a model, in which N real-scalar fields constitute a vector 



representation 



,<j) n ) of O(N). The case of N = 4 is of practical interest. Specifically, through computing 



the effective action, we analyze how the phase transition proceeds. The Lagrangian of the theory reads 



(d$ B ) 2 - fi 2 J 2 ] - ^{$1? + H B Ce- $ B ) 



(6.1) 



where fJL 2 B < 0, so that the theory describes the system whose ground state is in a broken phase in the classical limit. 
In ( |6.ip , we have introduced the term _ffs(*e • 4>b), with e an unit vector, which explicitly breaks O(N) symmetry. 
Noticing the fact that an UV-renormalization scheme for the symmetric phase (fi 2 B > 0) works |^7j as it is for the 
symmetry-broken phase (fi 2 B < 0), we introduce UV-renormalized quantities, 



\fZ(f>, n 2 B = Zf.fi 2 



ZxX. 



In terms of renormalized quantities, £ reads 



= \ [{d$f ~ ^V 2 ] - ^($ 2 ) 2 + HCe- $) + \{Z 1) [{d$) 2 - fiH 2 
-(Z x Z 2 -l)±($ 2 ) 2 +Afi\ 



(6.2) 
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where H = \J~ZHb- In (6.2), according to [ 
to make the effective action finite. 

For the purpose of later use, we here construct 



we have introduced the c- number counter term A/j, 4 , which is necessary 



£{(j>{x);0{x)) = £{^{x)+(p{x)) 



dC{0(x)) 
d(f(x) 



(6.3) 



To avoid too many notations, for C((j>(x);ip(x)), we have used the same letter "£" as in ( p.l\j and (|6.2|). As will be 
discussed below, <p{x) in (3.3) is the (classical) condensate or order-parameter fields and <fi(x) is the quantum fields, 
which describes the fluctuation around <p(x). Straightforward manipulation yields 



£((f>(x);ip(x)) 

^cond 



£' 



cond 



£' 



a 



£cond 

l 

2 
1 

2 
1 

2 
1 
2 

+ 3! (V ' 



£' 



cond 



mV 2 



(Z-l)(d0) 2 

(d$f - nS 



a + a 



1 



A 



{Z^Z — l)fl 2 (fi 2 — (Z\Z 2 — l)^y((/? 2 ) 2 



x_ 

~ 12 
1 
2 
A 
4! 



+2(v-0) 2 l 



{Z ^Z - 



- (Z X Z 2 



3! 
1) 



V, 



(6.4) 
(6.5) 



2 f 



12 



{^V 2 + 2(V-<?) 2 } 



2\ 2 



(6.6) 



It should be noted that £ enjoys 0(N — 1) symmetry. 

For obtaining an efficient or rather physically-sensible (perturbative) scheme, as in Section VE, it is necessary 
-|30|) to introduce weakly ^-dependent masses, M 7T (x) and M a (x), for which we assume that 

\^-M 2 (x)\, \M$(z)-MZ(x)\ = 0(S). 



(6.7) 



How to determine M 2 (x) and M 2 (x) will be discussed in subsection D. We then rewrite £((f>(x); <p{x)) in the form, 



Here £ c 



£((f>(x);tp(x)) = £ CO nd + £'cond + ^0 (<f>(x) ', 0{x)) + £ int + £ r 

and £' cond are as in ( |6.4[ ) and (6.5), respectively, and 

1 t 



r 

t^rc 



l 4>(x) A~ 1 (x) 4>(x) 

x fit- ftp- A 



3! 



4! 



£' 



: l kx) [x*(x)P„{x) + X a{x)V a {x)] ${x) 



(6.8) 



(6.9) 



(6.10) 



Here £' rc is as in ( |6.6| ) and 

A" 1 ^) = - 



A 



8 2 + M 2 (x) + -0 2 {x) ) P,(x) -\d 2 + M 2 (x) + -0 2 {x) ) P CT (x), 



A 



X«(a;) = Mf(z)-/x 2 (£ = tt, o), 
P»(a;)=I-| ¥ >(a;))( ¥ >(x)|, 

P CT (z) = (6.11) 

In these equations, the boldface letters denote N x N matrices in the real vector space: (Pa) is the projection 
operator onto the n- (cr)-subspace; I is an unit matrix. \(f(x)} is a unit vector along <p(x), and (<p(x)\ is an adjoint of 
\tp(x)). As has been noted above, when compared to M 2 (x) (£ = tt, <r), Xs( x ) is °f higher order. 

As mentioned above, the quantum fields </> describe the fluctuation around the condensate fields <p(x). The con- 
struction of perturbation theory based on (6^8) starts with constructing the Fock space of the quanta described by <j>, 
which is defined "on <p(x)." As stated in Section I, dealing with the dynamical aspect of the system, we are assuming 
that microscopic or elementary reactions can be regarded as taking place in a single spacetime cell. While, in a single 
cell, relaxation phenomena are negligible. This means that the theory to be developed in the sequel is consistent only 
when 0{x) does not change appreciably in a single spacetime cell. We shall assume that this is the case and treat x 
of (p{x) as macroscopic spacetime coordinates. 
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B. Effective action 



Spacetime evolution of the condensate fields <p(x) is governed by the Euler-Lagrange equation, which is derived from 
the CTP effective action [ip\. If one restricts to the case where <p(x) is space-time independent, ip(x) — > 0, Tc[<p\ 
reduces to minus the CTP effective potential, rc[</?] = — Vc (</?)■ The physical part of the CTP effective potential (cf. 
below) Vphysi'f) is the minimum free energy density of the system under the constraint, (4>(x)) = tp, the minimum free 
energy density with quantum and medium effects taken into account. A computational algorithm of effective potential 
in vacuum theory is established in ]3l[ | . On the basis of a generalization of this algorithm to the case of equilibrium 
thermal field theory, a calculational scheme of the effective potential and its physical part have been settled (§2] (see 
also H). 

We generalize this algorithm to the present case of out of equilibrium systems. Following standard procedure g] of 
obtaining an expression for the CT P effective action, we start with the CTP generating functional for Green functions 
in the interaction picture (cf. ( |3.7| )), 



Z C [J 



= JW C [J] _ 



(Tcexpi / C\${x)) + J(x) ■ ${x) ) 



(6.12) 



where, as in Section II, C (= C\ © C<i) is the closed time-path in the complex ti me-p lane a nd £ stands for the 
perturbative part of the Lagrangian density, which is to be identified below in ( 6.20| ). In ( |6.12 ), Wc[.7] is the 
generating functional for connected Green functions. Expectation value of </>{x), {<f>{x))j (= <p(x)), in the presence of 
the external source J(x) is obtained through 



tp(x) = 



SW C [J] 



8J(x) 

The CTP effective action Tc is defined as the functional Legendre transform, 

r c [<p\ = w G [J)- [ J{x)-(p{x), 

Jc 

from which the Euler-Lagrange equation follows: 

5T c r 



5ip(x) 



= -J(x). 



(6.13) 



(6.14) 



(6.15) 



The CTP effective action Tc is the g enera ting f uncti onal for connected, one-particle irreducible, amputated Green 
functions. Substitution of ( 6.12 ) and ( 6.15 ) into ( 6.14 ) yields 



(T c exp 



£{${x)) 



STr 



6<p(x) 



{<j>(x) - <p(x)} 



Redefining <p(x) — <p(x) as new quantum fields <p(x), we have 



(T c exp 



£($(x) + 0(x)) - 



5Tr 



8(p(x) 



fa) 



(6.16) 



Thus, the new quantum fields <f>(x) describe the fluctuation around 0(x). Here, for the purpose stated below, we have 
inserted the factor i5 (=1). Going to the two-component representation, we have 



ir[<p lt <p 2 ] _ 



= (r c exp 



1 / d i xJ2\(-) j - 1 ^ , (^)+^))- 



j'=i 



6f\<£l, 02\ 

80 j (x) 



<j}j(x) 



(6.17) 



where <pj and tpj (J = 1, 2) are the type-j fields and 



6<pj(x) 6<p(x) 



(i = M). 



(6.18) 
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The minus sign on the RHS of ( 3.1§D with j — 2 comes from the fact that, on C2, the time ticks away backward from 
+00 to -00 ( cf. 

Equ ation ( |6.17| ) is of the form which allows for computing f [^1,(^2] perturbatively. Perturbative calculation of 
( |6.17[ ) goes by taking suitable 



2 2 

3=1 3=1 



(6.19) 



being bilinear in quantum fields <pj (x), as the free hat-Lagrangian density, which defines 2x2 propagators of the fields 

<t>( x) = (0i (x ), 02 ( ^))- The form of C{<f)(x) + <p(x)) is already been obtained in the last subsection and Cq is defined 
in ( |6.9| ) with (6.11). The 2x2 matrix propagators will be constructed in the next subsection. The vertex factors can 
be reads off from 



<c-EH i-1 A0iC«O)-4) 
3=1 



\- st[0 1 ,0 2 ] 



3=1 



(6.20) 



Note that, to zeroth order (tree level), 

f(°)[^ 2 ] = ^2n j - 1 j d^xd^x)). 



(6.21) 



3=1 



Then, the exponent in ( |6.17j ) with for f on the RHS is free from the term being linear in 0, which means that 
the "one-point vertex" is absent. [In constructi ng C( 4>(x); <p(x)) above, Eq. (3.3), this has already been taken into 
account.] Thus, in terms of Feynman diagrams, ( |6.17 ) consists of bubble or no-leg diagrams. The loop expansion is 
defined as an expansion by S ((5-expansion) , which is set to 1 at the end of calculation (see, e.g., |34j]). Computation of 
r should be carried out successively: In computing Zth loop-order contribution, the contributions to T up to (I — l)th 
loop-order should be substituted for f on the RHS of ( 6.17 ). This part, — iS^ 1 Y^j=i $j ' ^/Stpj, plays a roll of 
eliminating the contribution of one-particle reducible diagrams. 

The one-loop contribution to T will be dealt with in Section VII. Computation of a contribution r( n ) from a multi- 
loop diagram, which includes n (> 1) vertices, goes in standard manner. Each vertex consists of two types, the type-1 
and type-2. Then, r( n ) is obtained by summing over types of vertices, so that T^ n > consists of 2™ terms, which we can 
write, with obvious notation, as 



p(n) 



f 1,^21 



E 



i 1 ,i 2 ,-,in [<Pl><P2\ 



(6.22) 



where ij (J = 1, 2, n) denotes the type of vertex. 

In the two-component representation, the Euler-Lagrange equation (|6.15|) reads (cf. above after ( 6.18 )), 



i-y-'jjix) c? = 1.2). 



In the limit of absence of the external sources, J\{x) — J2(x) — 0, 0\{x) — (p%{x) holds and ( [5.23 ) becomes 



5<pj(x) 



(7 = 1,2). 



(6.23) 



(6.24) 



The equation with j = 1 is the equation of motion for the physical or type-1 condensate fields, while the one with 
j = 2 is for the unphysical or the type-2 condensate fields. The physical effective action Y p h ys and the unphysical 
effective action T unp f lys are defined through 

ST phys [ip\ _ ST[0i,02. 



80{x) 

^1" unphys [' 
80{x) 



difl(x) 
6T[0 1 ,02_ 

Sip 2 {x) 



(6.25) 
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With suitable (functional) integration constant, ( |6.25| ) determines T p h ys and T unp i iys . In the case of thermally equi- 
librium system, the physical effective potential (T p hys with constant <p) has been computed in |]32"f up to two-loop 
order. 

At the previous subsection, in defining Co, we have introduced (arbitrary) "mass functions" M^(x) (£ = ir, a). The 
free hat-Lagrangian density in the present CTP formalism is given by (|6l9| ). Then, as in the fields, there are two types 
of Mf{x), the type-1 and type-2, which we write as M^(x) (j = 1, 2). [M^(x) = M|(x) with x £ C j (j = 1, 2).] As 
will be discussed in subsection D, various methods of determining Mh(x) are available. All those methods are based 
on the notion of "optimization" of some physical quantity. The determined M?j(x) generally depends on ipi(x) and 

0i{x). Let 02] Mj^i) M% 2 , M%i, M% 2 ] be the effecti ve ac tion computed to some order of loop expansion. Here a 

question arises: In computing the equation of motion (6.24), does the functional derivative 5/6<pj(x) (j = 1,2) act 
on Mf^ix) and M| 2 (x)?. The answer is as follows: If we intend to "optimize" the effective action, from which the 
equation of motion is derived, 6/6tpj(x) should act on M| 1 (x) and Al^ 2 (x). On the contrary, if we intend to "optimize" 
the equation of motion, 6/6<pj(x) should not act on M| x (x) and M? 2 (x). We take the stand that the effective action 
is an intermediate device to derive an equation of motion for <fj(x) (j = 1, 2), etc., and then the equation of motion 
is of physical importance over the effective action. On the basis of this observation, following |lQ| , we take the latter 
throughout in the sequel, i.e., 5/8(pj{x) does not act on M| 1 (x) and M| 2 (x). 



C. Propagator 

1. Retarded/ advanced propagator 



Here, we obtain a retarded/advanced propagator, A^mg, which is an inverse of 

A-\x)=A- 1 (x)t 3 , 



with A 1 (x) as in ( 6.11 ), under the retarded/advanced boundary condition (cf. fl4.2|) and (4.4)) 
We start with finding an inverse of A -1 (a;) under the retarded boundary condition: 



A-\x)A R (x,y) =l5\x~y). 



To this end, we introduce 

.(0)/ 



A%(x,y) = A%\x,y) + A%'(x,y), 



-iP-(x-y) _ 



P 2 -M\(X) 



where X = (x + y)/2 and 



A 



A 



MUX) = M*{X) + - 0\X), MUX) = M 2 {X) + - 0\X) 



(6.26) 



3.27) 



(6.28) 



(6.29) 



Mi{X) in (|6.28D should be understood to be M\(X) - ie{p )rj. 
Straightforward manipulation using (C2) in Appendix C yields 

74; 



A- 1 (x)Af(x ) y) = I6 4 (x-y)-2i J ±JL e -^^P^ 



dXv 



P 2 -Ml(X) P 2 -Ml(X) 
With ( |S3C|) in hand, we are now ready to obtain Ar(x, y) in (|6.27|) to the gradient approximation: 



(6.30) 



A R (x,y)^ J2 P d*) J 



-iP-(x-y)_ 



7T, tx 



P 2 -M\(X) 



-iP-{x-y) 



p-d x p a (x) 

P 2 -M 2 (X) 



P 2 -Ml(X) P 2 -M 2 a {X) 
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Using ( |C2D in Appendix C and Fourier transforming the resultant expression, we obtain 

p 2 

1 1 



A «^ P ^ £ P^M|(X) 



-iP 



(V(X)) 5^ (^(X)| 



(6.31) 



_P 2 -Af 2 (X) P 2 -A4 2 (X) 
It can be shown that Aji(x,y), inverse Fourier transform of A.r(X; P), satisfies 

A R (x,y) A 11 (y)=I6 4 (x-y). 

Form for Aa(x, y) may be obtained in a similar manner. Thus we obtain for Adi a g(X;P), the counterpart to the 
Fourier transform of Adi a g(u — v) in ( |4.31| ), 



A diag (X; P) ~ A<°> (X- P) + A^(X; P) 



(i) 



3.32) 



A^L^P) = diag(AHX;P), -A^(X;P)j , 

= ]T P c (X)di^A«>(X;P),-A«>(X;P) 

7T, CT 

ee ]T P ? (X)A£ g (X;P), 

£=7T, CT 

ASi 9 (X;P) = diag(A«(X;P),-AW(X;P) 



(0), 



.(0)/ 



where 



A« A) (X;P) = -zP M (|^(X)) 



(6.33) 
(6.34) 



P 2 - MUX) ± ie(po)^ P 2 - X 2 (X) ± ie(p )r? 



(6.35) 



Adiag{X\ P) corresponds to A^^P), Eq. fl4.6| ), in the complex-scalar field theory. 
The following form for Adiag(x,y) is also useful for later analysis: 

A dl ag(x,y)= J2 P sW A £ fl (^y) P s(y) + A £ ff (^y): 

£=7T, CT 

— ^-^-^A^ (X;P), 



2tt) 4 
d 4 P 



3 -iP-(x- 



cdiag (AW(I; P)a£> {X ; P), -A« (X; P)A^(X; P) 



(6.36) 



Let us find the property of the operators |^(X))3x fl (</'(X)| and \ip(X)) dx„ (</?(X)|. As defined above, |</?(X)) is the 

unit vector along (p(X) in the O(X) vector space, so that ((p(X)\tp(X)) = 1. Then (ip(X)\d x ^\y{X)) = (y(X)| d X „ 
\<p(X)) = 0. From these relations, we obtain 

P ? (X)(KX))dx>(X)|)P c (X) : 
P S (X) (\if(X)}dx. (<p(X)\) P C (X) : 
Thus, A d ^ ag in (6.32) and A^ s in ( |6.36 ) represent mixing between the 7r sector and the u sector. 



|^(X))dx>(X)| ife = aandC = 7r, 
otherwise, 

\<p(X))d x J<p(X)\ ife = 7ra n dC = a, 
otherwise. 



(6.37) 
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2. 2 x 2 matrix propagator 
Let us introduce "Bogoliubov" transformation (cf. ( |4.15| )), 



d 3 u ^ Bi r) (x, u; x )<f RA (u, x ) 



r=± 



0(y)= / d 3 «^*4 A (v, 2/0 )B^ ) (v,y;y ). 



As seen above after (5.6), £ g enjoys 0(N — 1) symmetry. We assume that the density matrix enjoys 0(N — 1) 

- (t) 

symmetry in the sense that the N x N matrices take the form, 

^i(Ji)( x ' u ; x o) = P«(x,2:o)Bi T ( ) i)' C) (x,u;a; )P c (u,a; ). 

4, C=T, CT 



The condition that the canonical commutation relations are preserved (cf. (4.13)) yields a solution, 



B 



, iR) (x,u;a;o) = P e (x, x )i?[ T ( ) ^ ) (x, u; x )P 4 (u, x ), 



where -B^"^' 1 is as in (4.16) with /( T ) — * • Translated into the |po| prescription (cf. ( [4.3l| )), we obtain for the 
propagator, 

A(x,y) = [ d i ud i vB L (x,u)A diag (u,v)B R (v,y) 



Bf{x,u) 



d^ud^v Y, P d^Hx,u)P i (u)A diag (u,v)P c (v)B^(v,y)P c (y), 
S 4 (x-u) fd x ' u ) 



5 A {x~u) 5 4 (x~u) + ft(x,u) J ' 
6 i {v-y) 5\v~y) 



(6.38) 



Substituting (6.36) into (6.38), we obtain 

A(x. y) — Jd\d 4 v Y Pd^\x 1 u)P ( (u)A^l g (u,v)P ( (v)Bf(v,y)P^y) 

^— 7r, a 

+ fd\d i v ]T P«(^ ) (^^)P?(")A^ g K«)P C («)4 C) («,2/)Pc(2/)- 
From the first term on the RHS, we pick out 

p s Oz)PeWP £ («)Pe(y)- 

Using P^(u) ~ P^(a;) + (u — x)^dx l ,P^{x) 1 etc., we obtain 

Eq. O = P«(x)P ? (y) + (« - a:)pP e (x) (c^P^)) P ? (x) 
+(v-y) fi P e (y) (d y Ps(y))Pt(y). 

In the second (third) term on the RHS, we have set P^(w) — P{(y) — Pf(z ) (P^(a;) ~ P^(u) — P^(y)), which is 
allowed since the term already contains the .X^-derivative. The relations ( 6.37 ) tell us that the second and third terms 
on the RHS vanish. Thus, we obtain, for the first term on the RHS of ( 6.39| ), 



(6.39) 



(6.40) 



]T P c (aO J d i ud\B^\x 1 u)k^l g {u 1 v)Bf{v 1 y) 



(6.41) 



35 



The quantity in th e square brackets is the sam e as i n comp lex-sca lar the ory dealt with in Section IV. Fourier trans- 
formation of ( 6.41 ) may be ca rried out using (|6~38l) , ( |oq ) and (|c|) - @. The last term in may also be 
computed by using ( |6.36j ) and ( |6 .37 ). Fourier transforming the resultant expression, we have 



+2iP>(X)) 



a r^r 
AM - A\A A -A A A A 



{f a A R A R - UA\A° A )A, 



d{ V {X)\ 



where A R = A R (X;P), f a = f a {X), etc. 

After all this, we finally obtain for the Fourier transform of A(x, y) on x — y: 

A(X;P) = A (0) (X;P) + A (1) (X;P), 
A(°)(X;P)= J2 P dX)B^(X;P)A2 g (X;P)B ( i\X;P) 

£= 7T, 0" 

= J] P € (X)A«)(X;P), 



where 



Here 



£=7T, 
1 



Pi f) (X;P) = 
A«)(X; P) = A^(X; P) + A^' (X; P), 



1 1 + f^X-P) I " Bf{X;P) - 



_(l + f s (X;P) fs(X;P) 
1 1 



A^i(X;P) 



>(€) 



,(f) 



a(«) 

A«>' (X; P) = / £ (X; P) [A«> (X; P) - A«> (X; P)] i+, 

AW(X;P) = -1 ^ P ? (X){/ e (X;P),P 2 -M|(X)}[(A^) 2 + (A« ) ) 2 



-4 4 



dP a 



(Ag } 



f2 = f) +(/ CT [(A^) 2 -2A^A i? 



(A^) 2 ] - /.[(iWf - 2AWAW - (a^) 2 ; 



n = 



(A«-A^) 2 

(aW-a('') j -(aW-aW) j 



-(A^-A^) 2 



Note that is related to the "genuine" number density through (cf. ( [4.3C| )) 

ft(X;P) = 0( Pn )N^(X;E^;p) - (1 + 9(- Po ))N® (X; E®; -p) 



(6.42) 
(6.43) 

(6.44) 

(6.45) 

(6.46) 
(6.47) 



(6.48) 



(6.49) 
(6.50) 



3. Counter Lagrangian and renormalized scheme 



Following the procedure in Section VE (cf. (5.54) to ( |5.55| )), we can extract the counter hat-Lagrangian density C c 

through analyzing A _1 (x) A(x, y) and A(x, y) A -1 (y). £ c thus found gives a two-point vertex function V C (X;P). 
Straightforward manipulation yields 
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= J] P«W{/«(X;P),P 2 -VW|(x)}i_ 



(6.51) 



(6.52) 



Let us adopt the physical-/ schem e. In the case of complex-scalar field, we have adopted the "renormalization 
condition" V C (X:P) = (cf. ( 5.55 ) and ( 5.56| )). In the present case, we cannot set *V C (X;P) = 0. Then, we adopt 
the condition, 



{f c (X;P), P 2 ~Ml(X)}=Q, 



(6.53) 



which is the same as in the complex-scalar- field case. Then, in V c in ( 3.52 ), the second term on the RHS survives. 
Let us compute the V c -resummed propagator A c _ res . Since V c contains the A^-derivative, one can directly write 
down the Fourier-transformed form: 



A(X;P) 



l+^{-V c (X;P)A(X;P)} 



(6.54) 



In 2 x 2 matrix space, V c is proportional to A_, Eq. ( 3.16| ), which shows that, in ( |6.54 ), the contributions with n > 2 
vanish. Thus, after some manipulations, we obtain, 



*c—resum 



(X;P)~-2i(f 7T -f rT )P fM 



\<p(X)) dx. {v{X)\/Xl^ A + \y{X))d x M X )\^RXA 



A, 



Summing this to A(X;P) in (6.42), we obtain ( 6.42| ) - (6.48) with ( 6.53| ), provided that Cl in ( |6.4£ ) is changed to 



R A il^R 



= no+(/„[AS'-Alf'][A£' + A A 



-/„ [A%> - A£W + - 2A<f>] ) A + 



(6.55) 



This completes the deduction of bare propagators i n th e bar e-/ s cheme. The physical-/ scheme may be constructed 
in a similar manner as in Section IVC. In place of ( |5.27 ) or ( [5.29| ), we have, with obvious notation, 



X { $ f (X;P) = (£ = *,*), 
which lead to the generalized Boltzmann equations. 



(6.56) 



D. Determining M%(X) and M%(X) 



In this subsection, we present a scheme for determining so far arbitrary M|(A) (£ = n,a), introduced in (6.7) 
(3.11). Various methods of fixing M|(X) are available. Among those, we mention the following three methods 



1) Resummation method. Argument in Section V shows that the quantities to be resummed are the real part of the 
retarded and advanced self-energy parts, 

£ fl (X; R) = S„(I; R) + £ 12 (X; R) = (£ A (X; R))* . 

Determine M|(X) (£ = tt, a) such that the real part of X#(X;i?) computed up to required Lth loop order, 
Re ~Eft(X; R), vanishes at i? M = Rg with Rq a suitable four momentum: 

Ret^\X;R = R ) = 0. (6.57) 

2) The principles of minimal sensitivity ]35| ]. Total Lagrangian density C in ( |6.§| ) is independent of arbitrary 
functions M%(X) and M%(X). Let S be a quantity, in which we are interested. Then, if we could evaluate S to 
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all orders of perturbation series or loop expansion, it is independent of M?(X) (£ = n, a). On the other hand, S 



computed on the basis of £ in (6.S) up to the Lth order, S^ L \ depends on M|(X). This observation leads one to 
impose the condition 



6E^ 
SMl(X) 

r2i 



= 



(£ = 71", a), 



which serves as a determining equation for M^(X). Note that, for S = r p ^ ys or 5 = ST p h ys /S0(X), this method 
works for L > 2. 

3) XTie criterion of the fastest apparent convergence [ [36| . Compute 3 up to Lth order, (L > 2). Then, 

determine all the parameters of the theory, such as the coupling constant, mass functions M^(X) (£ = 7r, er), etc., so 
that all the terms but the first term of the series vanish. 

A number of works, which employ the methods 2) and 3) above for analyzing the equilibrium effective potential, is 
cited in 0. 

In this paper, we follow the method 1). The requirement in 1) gives the gap equation. 



E. Summary of this section 



In this section, we have found that the spacetime evolution of the condensate fields <p(X) is governed by the following 
system of coupled equations: 



El) Generalized relativistic kinetic or Boltzmann equation (6.56) 



E2) Gap equation (6.57). 

E3) Equation of motion for tp(X), (|J|) with j = 1 . 

As seen in Section VC, ( |6.53j ) in El) is a determining equation for the number-density function f^(X; P) (£ = tt, a) 
under a given initial data f^{X^ = Ti„,X;P). This equation involves yet unknown functions M|(A) (£ = 7t,ct). 
Equation ( |6.57| ) determines M|(X) in terms of f^(X; P). Thus, ( 6.53| ) and ( |6.57 ) interplay to determine M|(X) and 
f^(X;P). Through computing the effective action T using M|(A) and f^(X;P) thus determined, the equation of 
motion (6.25) for ip(X) is obtained. Solving it under a given initial data determines the spacetime evolution of <fi{X), 
which describes the progress of phase transition. 



VII. LOWEST-ORDER EFFECTIVE ACTION IN THE a MODEL 



In this section, we co mpu te the CTP effective action to the lowest nontrivial order, which comes^ from an one-loop 
diagram with £q as in (|6.9[) . The contribution is 



f = ^Trln^f 3 |l+ £ P.Mld-^ 



(7.1) 



where 'Tr' acts on matrices in spacetime, 2x2 matrices in the "complex-time plane," and N x N matrices in the 
O(N) vector space. In (O), 



(f 3 d- 2 ) (x,y) = -A(x,y) 



m = 



with A(x,y) as in (4.31) and 



3 The difference between f in (7.1) and f computed using the physical-/ scheme, whose form has not been displayed explicitly, 
is of higher order. 
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M\{x) = diag [M 2 (1 {x), -M 2 ^(x)] (£ = it, a) 

A 2 

(p = diag((^ 1 2 (a;), -<p%(x)) , 
Pt(x) =diag(P a (a ; ), -Pc 2 (a;)). 

Here ^ (i = 1,2) are the type-« condensate fields, .M 2 ^ and (« = 1,2) are as in ( |6.29 ) and (6.11), respectively, 
with (fii in place of (p. 

For extracting the physical effective action from (7.1), we use a method |32j of computing physical effective potential, 
in which the tadpole method is employed . The method starts with taking the derivatives, 



ST 
SMl(x) 



|trP ei (x) [ld' 2 (x,x) - {cT 2 (P«Ml +PaM 2 a ) d~ 2 } (x,x) + .. 

(£ = 7T,cr), 



(7.2) 



where "tr" acts on 2 x 2 matrices and N x N matrices. Then, we set 2 (x) — 9?| (x) = ip 2 (x). Equation ( |7.2| ) now 
becomes (cf. ( |6.25 ) ) 



ST 



SM 2 £1 (x) 



Vl = V2 = <fi 



= -r tTP s (x)G n (x,x) 



ST 



phys 



SMUx) 



(7.3) 



where 



G = -f 3 



Id 2 + P^Ml + P a Ml 



(7.4) 



V>1 = V2 = <P 



The one-loop contribution to the physical effective action T p h ys is obtained |32| by integrating (7.3) over Ai 2 (x) 
Incidentally, the "unphysical part" of T, which comes from 5T/5A4^ 2 (X), becomes 



unphys 



phys 



(7.5) 



So far, we have made no assumption on the x-dependence of Ai 2 (x) and tp(x). Here we recall that M^(x) and <p(x) 
weakly depend on x. We first note that 

Sp[dx^dX)} = Sppe^PcW] - (£,C = 7r,«7), 

where 'Sp' stands for the trace operation that acts on an W x JV matrix. Then, the term with dx^P^(X) does not 
appear. This means that, within the gradient approximation, we can assign common spacetime coordinates for all P^ 
in (7.4). Thus, including the zeroth-order contribution (cf. (6.21)), we obtain 



phys 



d i X j phys (X), 



lvhys{X) = \ [{d0{X)f - ^ 2 (X)f + H(*g. 0{X)) 

+\[{N-\y h (x) + la {x)], 

1S(X) = - [ dM\(X)mX) (£ = n,a). 



Here 



UX) = i J ^ i ^ i {X;P)=lf\x)+lf\x), 



if\x) 



lf{X) = 



d 4 P 



(2tt) 4 P 2 -M\+ir)' 
d 3 p N^(X;SP(X),p) 



(2tt) 3 



4°(x) 



(7.6) 

(7.7) 
(7.8) 

(7.9) 
(7.10) 

(7.11) 
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where £p (X) = yp 2 + A4^(X). Thus the physical effective action T p h ys is given by an integral of the "local effective 

action" JphysiX) over the macroscopic spacetime coordinates X^. Computation of "f p hy S (X) goes essentially in the 
the same manner as that of the physical effective potential. 

In the following two subsections, we separately analyze the positive-mass-square region and negative-mass-square 
region. As mentioned in Section I, our scheme does not apply to the region near the transition point between the two 
regions. 



A. The region M%{X), Ml{X) > 0. 



The lowest-order contribution to 1^(X), Eq. (7.9), has already been computed in Section VD (see (5.34) - (5.37)): 

M 2 AX) 



jf\x) = -- 



(4tt) 2 



2 , M\{X) 
— — In — ^5 — 



1 POO 

lf\x) = — / d£j&-M\{X) (N^(X;£)), 



where 



with dflp a solid-angle element in a p-space. Substituting these into (7.8), we obtain 



^(X)=^°\x)+^'(X) 



(P), 



with 



7 f(X)=i 



d 4 P 



In (P 2 - M 2 {X) + i V ) 



1 



7f(X) 



32tt 

1 

3^2 



2 M\{X) 



2 3 , M\{X) 



- + - - In 
e 2 



M 2 d 



M f (X) 



d£ (£ 2 - M 2 AX)fl 2 (N^(X;£)). 



(7.12) 

(7.13) 

(7.14) 

(7.15) 
(7.16) 
(7.17) 



2/e in ( 7.16] ) is canceled by the cou nter Lagrangian C' cond , Eq. (6.5). We note that, to the present one-loop order, we 
can replace — (Z^Z — l)/i 2 (p 2 /2 in (S.5) as 



(Z^-l)/xV72 



2(N + 2) 



(Z^Z-l) [(N ~ 1)M 2 + 3M 2 ] ^ 



(7.18) 



The difference b etwee n the LHS and the RHS, which is proportional to the UV-diverging factor 1/e, is of higher order 
(cf. above after ( 6.11 )). This difference is canceled by the higher-order one-loop diagrams, one of which includes the 
vertex coming from *0x 7r P 7r </)/2 (s C rc ) and the other one includes the vertex coming from t (j}\ a 'P a 4>/2 (s L rc ). Since 
we are choosing the MS scheme, we have 



Z = 1 



p 487r 2 e 



Zx = 1 + A 



7V + 81 
48tt 2 i' 



(7.19) 



In a similar manner, we obtain for A in (6.5), 



N 1 



32tt 2 e' 



After all this, we obtain for the local effective action, to one-loop order, 

j phys (x) = 7 ^ s (x) + 7 !f h > ys (x), 



(7.20) 
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where 



Ml(X) 



(8tt) 2 



(iV-l)A<(X)ln 



A 

4! 1 

MUX) 



M 4 a (X)ln 



4^) = , (iV-l) 7 W(X)+ 7 W(X) 



(7.21) 
(7.22) 



It should be emphasized that the factor of the u tp 2 term" in ( 7.21 ) is neither M 2 T (X) nor .M 2 (X) but /i 2 . 

In the case of thermal equilibrium system with constant ip, —"f p hy S reduces to the effective potential (density) of 
the system, V(ip) = —j p hys(<p)- Besides the terms due to the different UV-renormalization scheme, V(ip) obtained 
above agrees with the result of flOl , where M 2 (X) = M 2 (X), being constant, is assumed. 

Now we are ready to write down the equation of motion for the condensate fields (cf. Section VIB), 



50{X) 

n(0 2 (xj) 



= - [d 2 + Q{0 2 (X))] <p(X) +He = 0, 

= r, dlphys{X) 

d0 2 {x) 

= i? + ^0 2 {X) + ^[(N- l)l n (X) + 31 a (X)] 

= n o (<p 2 (x)) + np(0 2 (x)), 



(7.23) 

(7.24) 

(7.25) 
(7.26) 



where 



n o (0 2 (x)) 



np(0 2 (x)) 



» 2 + ^ 2 (x) 



96tt 2 



(N-l)Ml(X) In 



Ml(X) 



■3Ml(X) In 



MUX) 



(N ~1)1^{X) + 31^>(X) 



(7.27) 
(7.28) 



Here jf\x) (£ = tt, a) is as in (^T|). 

It is wo rth mentioning in passing that t he eq uation of motion for unphysical or type-2 condensate fields arc obtained 
from ( |7.5| ), which is of the same form as ( 7.23| ). Since the time argument of the type-2 field is on C2, the type-2 field 
"lives" in the "time-reversed world." Then, one can say that the equation of motion is time-reversal invariant. 

Since we are assuming that 0{X) does not appreciably change in a single spacetime cell (cf. Section I), (7.22) 
applies only in the region of <p(X) where |£1(</5 2 (X))| < 1/(L M ) 2 (/i = 0, 1, 2,3) with the size of the spacetime cell. 

Let us study the small oscillation of <p(X) at the vicinity of the quasiminimum point <f(X) — tpo(X)e: 



H. 



Setting 0(X) = ip Q (X)e + e(X), we obtain, 

d 2 + n^(^ 2 (x)) 

where e x (X) = e{X) - (*e • e{X)) e and 



H 



MX) 



( t e-e(X)) + d 2 MX) 
, h 



MX) 



MX) 



0, 
0, 



(7.29) 

(7.30) 
(7.31) 



n^(tf(x)) = n ( °\0 2 (x)) + rt;\^(x)) 



with 



d;\$(x)) 



\Mx) 



AW) 
18 



96tt 2 



(N-l) 



(N- 



l)ln^ + 9 1n^l 



0(X)= vo (X)i 



di n {x) , n di a (x) 



dMl{X) dM 2 a (X) 



(7.32) 
(7.33) 



41 



e±(X) describes quasi Goldstone mode. Equation (7.3C) may be used only in the region of <fio(X) where 
|^)(^ 2 (X))|<l/(^) 2 . 



1. "High-temperature" expansion 



The system is characterized by a sca le parameter or several scale parameters of mass dimension. Let us assume 
that, for small £, (N®(X;E)) in fl7.13|) behaves as 



(N^(X;£)) 



£ 



■C ( (X) for T K {X) >>£, 



(7.34) 



where T^(X), which we call the temperature, is proportional to the scale parameter or to some combination of the 
scale parameters. In the case of locally thermal equilibrium system, the averaged number-density function reads 



(N®(X;£)) 



1 



e S/T(X) _ I ■ 

where T(X) is the local temperature of the system. Then, we have 



(7.35) 



T t (X)=T(X), Ct(X) = --. 
Let us analyze j p hys{X) in the high-temperature region. We start with rewriting ( [7.17] ) as 

T f»(X)= 7 <f'(X) +7 («(X) 



d£ 



£ 3 ~ -£M\{X) 



4 -b 5 (X)T*(X) - -ct(X)M$(X)7?(X), 



d£ 



9{£- M S (X)) (£ 2 - M 2 AX)fl 2 -£ 3 + %£M\[X) 



x(N^(X;£)). 



(7.36) 



(7.37) 



In the case of locally equilibrium system, b^(X) = c^(X) = 1. At "high temperature," 7^ dominates, |7^| << 7^- 



7^2'' may be estimated by substituting (7.34) into (7.37) 



1 { £\X)^-M\{X)T,{X). 



7 f(X) ~ ^(X)T*(X) - ^{X)M}{X)T 2 {X) + ±-Ml(X)T^(X) 



Putting all this together, we obtain 
The local effective action ^ p hy S (X), Eq. ( 7.2C| ), now becomes 

l ph ys{X)^ 1 ^ lys {X)+ 1 ^ ys {X) 1 



(7.38) 



where lfh ys (X) is as in ( |7.2l|) , and 



1 

90 



I V - I i { nX(X)T*(X) - ^(X)Ml(X)T 2 (X) 
- \ ;r/, T i \ r^ X ) - l j-c a {X)Ml{X)T 2 {X) 



12tt 



(N -l)Ml(X)T^(X)+Ml(X)T^(X) 
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In the case where Ml(X) « T„(X) and M/^(X) » T a (X), the "cr-terms" here and below are dropped 



High-temperature expansion of fl/s((p 2 (X)), Eq. (7.28), and Clg (0$ (X)), Eq. (7.33), reads, in respective order, 



(N-l){ c w T^(X) - - M V (X) %(X) 

7T 



-3 <^ c a T 2 (X) --M a {X)%{X) 

7T 



and 



'M„(X) M a (X) 



2. "Low-temperature" expansion 



Here we consider locally thermal equilibrium system, which is characterized by (7.35). In the low-temperature limit, 
T{X) « M\{X), lf ] {X) in d7TT7t ) may be computes as 



Tf>(*)~-ij / dE(E 2 -M\(X)) 
3?r 2 J M({ x) 



3/2 



oE/T{X) _ 1 



1 /*°° 

(Ml(X)) 3/i e -Ml(x)/Tlx) I dxx V2 e ~,/T { x) 



(M 2 (X)f 4 (T(X)f 2 e -M ( (x)/T(x)_ 

v 2ir 



Then, Jp^ vs (X) in (7.22) is approximated as 



L- (T(X)f 2 \(N 1) (Ml (X)f 4 e-MAxynx) 



2 3 / 2 7T 

+ (Ml(X)) 



3/4 & -M a {X)/T{X) 



Qf3((p§(X)) in (7.28) becomes 



Slp(0 2 (X)) 



(T(X)f 2 \(N - 1) {Ml(X)f A e -MAx)/T { x) 



and Ug\(pQ (X)) in ( [7.33 ) becomes 



n£% 2 p0) 



12V2tt 

+3 (Ml(X)) 1/A e- M ~M T W 



-^EMl ( T (X)f 2 \(N - 1) (Ml {X)Y 1/4 e -MAX)/T(x) 
72v27r L 



+9 (A^X))^-^*)™ 



(7.39) 



B. The region M|pf) and/or A-t^pf) < 



It can be seen that, through analytically continuing the one-loop effective action obtained above to the region where 
Ml(X) and/or Ml(X) < 0, an imaginary part emerges and then we face the issue of interpreting it. In the case 
of effective potential in vacuum theory, this issue has been studied, e.g., in fl37|| . Let us summarize the observation 
in p7| . The quantity analyzed there is essentially (7.15). In (7.15), the quantum modes with p 2 > \M 2 (X)\ remain 
harmonic oscillator as in the case of M^(X) > 0. However, the long-wavelength modes, p 2 < \M^(X)\, have imaginary 
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frequencies and thus unstable. These unstable modes turn out to develop imaginary part of r y^ t) \x) 1 Imj^ °\x), 
the emergence of which is a signal of appearance of perturbative instability. The 2m 7! (X) may be interpreted as 
half the decay rate per unit volume of a particular quantum state. More precisely, 2m 7^ (X) is the sum of half the 
decay rate per unit volume over all unstable modes, 4>iong{x), which are localized at (4>iong{xo = (xo)mi x )) — at the 
initial time Xo = (xo)m- Here 4>i on g{ x ) stands for the long-wavelength components of the quantum field 4>(x). The 
perturbative instability arises due to the growth of {4>iong{x)) ■ Owing to this instability, the original quasiuniform 
configuration beaks up into many uncorrelated domains. Within a domain the quantum field is correlated and 
{4>iong{x)} is increasing in time. From one domain to the next, (4>iong{x)}'s vary randomly from positive to negative. 
These domains is expected to have size at least as great as |.M 2 (X)| -1 / 2 , the minimum wavelength for an unstable 
mode. 

As seen above, leading one-loop contribution to the effective action is given by an integral of the "local effective 
action" "f p h ys (X) over the whole macroscopic spacetime region. We recall that 7 P / lys (X) is a function of X that 
labels the spacetime cells. Then, in a single cell, X of 7 P h, 2/s (X) is approximately constant. Moreover jphysjX) is 
the same in form as the effective potential (density). From this observation, together with the conclusion BlJ for the 
zero-temperature effective potential, we assume that the imaginary part of the effective action in the present case also 
reflects the proper physical situation. 



For the purpose of performing analytic continuation, we return back to (7.14) - (7.17) with A4 2 (X) > 0. Equations 



( 5.34 ) and ( |7.15| ) tell us that A4 2 (X) in ( 7.17 ) should be understood to be M^{X) — irj with rj a positive infinitesimal 
constant. Now let us analytically continue the function 7^ = 7^ + from the region A4 2 (X) > to the region 
A4 2 (X) < 0. The UV-renormalized vacuum part 7t , Eq. ( 7.16| ), goes to 



(0) 



32tt 



2 M\(X) In 



-M\{X) 



32tt 



M\(X) 



(m 2 jx) < 0). 



(7.40) 



It is worth mentioning that a "direct calculation" of ( [7.1C ) yields the correct result for Im 7^ 

d 4 P 



(0). 



iImlf ] {X) =ilm 



dp 



4tt 2 



^p 2 -M 2 (X) 



from which we obtain 



i Im 



7£ 



(0) 



J dMl{X)Hmlf\x) 



1 

32^ 



Mt(X), 



(7.41) 



in accord with ( |7.4C| ). 

As to 7^, we rewrite (7.17) as 



(/?) 



6tt 2 



d£(£ 2 -M 2 JX)+t V )(M 2 (X)-£ 2 - lv ) 1 / 2 (N^(X;£)) (M 2 JX) > 0), 



Ce 



(7.42) 



where (N®(X;£)) is as in ( |7.13[ ). In a complex-^ plane, the function (A4^(X) — £ 2 — i?/) 1 / 2 has branch points at 



£ = ±£b r with £br = (A4 2 (X) —ijf) 1 / 2 . We take a branch cut associated with +£b r to be a straight line [£br, +00- 



ta) 



where a = —Im(A4 2 (X) — irf) 1 ^ 2 . The integration contour in ( 7.42 ) is defined so as to enclose the branch cut 
[£br, +00 — ict) in a clock- wise direction, 

Cs : +00 - i(a + e) — * M$(X) - e - i(a + e) 
— > M^(X) - e + i(e - a) — > +00 + i(e - a), 
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where < e < a. 

Continuing to the negative A4%(X) region, the branch points ±Et, r go to 

± £br ± (-li-MliX)) 1 ' 2 + a) (a = 7m(-X|(X) + it?) 1 / 2 ) 



(7.43) 



and the branch cut associated with £ br goes to [£b r , (1 — i)cz] ® [(1 — i)a, +oo — ia). The integration contour Cs is 
deformed accordingly. We assume that (N^(X;£)) is analytic in the region 

a - e < Re£ < (-Al|(X)) 1/2 + a + e, 
-(-Mf pf)) 1 / 2 - e < Im£ <e-a. 

Then, we can deform the branch cut associated with £b r to 

[£ br ,£ br + {-MliX)) 1 ' 2 ] [8^ + (-Ml(X)) 1 / 2 , {-MUX)) 1 ' 2 - ia] 
8 [{-MUX)) 1 ' 2 - ia, +oo - ia). 

The integration contour Ce is deformed accordingly. 
Then, we obtain, after some manipulations, 



if = i /„. d * - Ml{X)f' 2 (N^(X; z)), 



(7.44) 



where C = C[ ® C 2 ® C 3 with C{ = -i{-M 2 (X)y/ 2 -> (1 - i)(-M|(X))V 2 , C£ = (1 - i)(--M|(X)) x /2 
{-M 2 AX))V 2 , and = (-.M 2 ^)) 1 / 2 _ +00 . 



After all this, we obtain, for the local effective action, 



l P hys(X)=J^ ya (X)+J^ y8 (X), 



(7.45) 



where 



^7&(^) = 5 [{d0(X)f - n 2 2 {X)} - £ (v 2 (X)) 2 + hCe- 0(X)) 



1 



(87T) 2 
1 

64tt 
1 



4! 

(TV - 1)M*(X) In '• M 2 ; 3 (X 2 )I + M*(X) In 



|Af 2 (X)| 



e 2 / 3 ^ 2 



^7i( X ) = ^ K N ~ l )° i-Ml(X)) Mi(X) + 6 (-Ml(X)) Mi(X)] 



(N-l)^\X) +1 (f)(X) 



(7.46) 

(7.47) 
(7.48) 



As mentioned at the beginning of this subsection, ImT p hy S {X) causes the perturbative instability, which in turn 
leads to formation of a domain structure. In most cases analyzed so far in the literature, imaginary part appears in 
the CTP effective action V in the form, 



N 



ImT = i^2 J d 4 xd 4 y(ipi(x)-ip2(x)) l "{ l: i(x 1 y)((pi(y)-ip2(y))j, 



(7.49) 



where 7^ are real functions. Recalling that the limit 02{x) — > <Pi(x) is taken at the final stage, we see that (7.4£) 
describes small fluctuations. In the present case, (7.5) shows that ImT under consideration is not of the form (7.4S) 
and describes "large" fluctuations.^ Here, for dealing with ImT p h ys , we follow a standard procedure ]3q |, which starts 
with writing exp iT p h ys in the form, 



Even in the present case, the small fluctuations of the type (7.4S ) also appears, which we do not consider in this paper 
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cxpi J d 4 Xj phys (X) 

i / d 4 X Re^p h y S (X) - I d 4 X Im-f phys {X) 



exp 



Provided that 



[dip] exp 



d 4 X Imy phys (X) 



< oo, 



(7.50) 



(7.51) 



we can make a functional Fourier transformation of exp[— J d 4 X Im^j p h ys {X)\. 

cxp[- / d 4 X Im^ phys {X)] = exp - / d 4 X Im-y phys (X) 

J l J = o. 

x /[dC]exp[- / d 4 XG{Q 2 {X)) + i I d 4 XX{X)-p{X)]. 



(7.52) 



Substituting (7.52) into (7.50), we obtain 



exp 



d 4 X Imj phys (X) 



ip = o. 



x(exp 

where the ^-average (...)^ is defined as 



i / d 



A X [Re lphys {X) + f(X) ■ 



(...) c = /[rfC]e-/ d4XG(C " (X)) 



(7.53) 



Equation ( [7.53[ ) tells us that Rej phys (X) + *<f(X) • ^(X) may be interpreted as the local effective action, from which 
follows the equation of motion for ip(X), 



dRe-f phys (X) 
d0(X) 



+ CW = o. 



(7.54) 



Thus, the auxiliary fields ((X) represe nt ran dom forces, whose character is determined by G(( 2 (X)). 

Taking the functional derivatives of ( 7.52 ) with respect to ip and then setting (p = 0, we obtain the relations: 



(i)c = i, 
(C(x)) c = o, 

(UX)C 3 (Y)) c = 26 t3 6 4 (X-Y) 
= 5 l3 S 4 (X Y) 



dIrrvy p h ys {X) 



d<p 2 (X) 



= 



A {(at _ 1)0 (-MUX)) MUX) 



+3 9(-M 2 a (X)) M 2 (X)} + (N-1)0{-M 2 (X)) dIm ^( X ) 



d<p 2 (x) 



0=0 



dip 2 {X) 



= 



(7.55) 



(Ci(X)Q(Y)<: k (U)Ci(V)) (c 

= ux)Cj(Y)a(u)Ci(v)) c - [<c ( (x)(,-(y)) c (a(u)Ci(v)) c 



4G 



= -2 <5 4 (A - (7) <5 4 (V^ - r)<5 4 (X - Y) [Sik Sij + 5 kj S u + S iS S lk ] 
A 2 



576tt 



[ {N _ i) (-Ml(X)) + 90 (-M CT 2 (X))] 



+ {N-l)6(-Ml{X)) 



d 2 Im^\x) 



d(^(X)f 



= o 



;-m ct 2 (x)) 



d 2 Im^\x) 



d{^{x)Y 



<P = o 



(7.56) 



etc. As mentioned in Section I, the macroscopic spaceti me c oordi nates are the indices that labels the spacetime cells, 
whose size is L» (,u = 0, 1, 2, 3). Then, S(X^ - Y^) in (|7.55|) and (|7.56|) should be understood to be 



8{X» - F M ) 



(no summation over /z), 



where etc., are the ^-component of the label for the spacetime cells. 

Let 0{X) = ip(X; [£]) be a solution to (7.54) under a given initial data (at the initial time Xq — Tm), 0{Xq — Ti n , X) 



and d<p{X ,X)/dX 



X = T in 



Then using ( 7.55] ), ( 7.56 ), etc., one can compute the £-averaged correlation functions 



etc. 

A multi-point correlation function consists of "connected" and "disconnected" parts. For example, a two-point 
correlation function (<p i (X)(pj(Y)) l ^ takes the form 

where the first term on the RHS stands for "connected" part and the second term stands for "disconnected" part. It 
is convenient to introduce a generating functional 



(7.57) 



which generates ((pi 1 (X 1 )...ip in (X n ))^ through 



(< Pil (xi)... iPin (x n )) < = Hrn 



l^Sy^Xj) 



Q[vM 



y = o 



(^-averaged "connected" correlation functions are obtained through 

5 



(^ 1 (X 1 )...^„(X n )) C c = (-i) n n 



8y i:j (Xj 



■lnG[y,<p\ 



v = o 



(7.58) 



It is worth mentioning here the equation of motion for type-2 or unphysical condensate fields. From (7.5), it is 
obvious that, for the unphysical condensate fields, all the above formulae apply with the substitutions, 



'Yphys * 'Yphys 



Imj phys — > -Irwyphys- 



This means that ( 7.55 ) and ( [7.56] ), etc., change sign. Recalling again that the type-2 field "lives" in the "time-reversed 
world," we see that the random forces that statistically upsetting the system act in completely opposite way in the 
time-reversed world. 
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1. "High-temperature" expansion 



In a similar m anner as i n th e case of Ai'j(X ), M% ( X) > 0, one can obtain a "high-temperature" expansion of 
jphys(X). Using $j\34|) in ( ^g ), we obtain (cf. ( |7^ ) - ( p8| )) 



1 



C" 



£ 3 - -£\^(X) 
(z 2 -.Mf(X)) 3/2 - 



rfz 



^b 6 (X)Tf(X) ±ct(X)Ml(X)T?(X) + ^ (-M 2 (X)) 3/2 T 4 (X). 



Gtt 



(7.59) 



It is to be noted that (7.59) is directly obtained by analytically continuing the asymptotic form (7.38) from the 
region M^(X) (> 0) to the region M^(X) ~ ir\ (M^(X) < 0). It is also worth mentioning that the following "direct 
calculation" reproduces ( 7.59| ). We start from ( |7.37| ), 



dp ■ 



o ^+MI(X) 
J™d£ (s 3 -^£Ml(X)\ (N®(X;£)) 



(N®(X;Jp> + M*(X)) 



which is valid for A4 2 - > 0. At this stage, let us naively change Ai^(X) from positive to negative. Then we have 



dp ■ 



jp>-\M\{X)\ 
J™d£ U* - \EM\{X)\ (N®(X;S)) 



(N^(X;^-\Ml(X)\) 



Substituting (7.34) and using ( 7.36| ), we obtain 







3tt 2 
r 00 


Jo 


J^/\Ml{X)\ 


iT^X) 


p CO 


3tt 


/ ° 

Jo 



dp ■ 



dpp\ x /p?-\Ml(X)\ + - 



dppU{p 2 -\M\{X)\) 



\Mj(X)\ 



^-\MI(X)\ : 



Gtt 



\Ml(X)\ 3 /%(X). 



This together with (7.36) reproduces (7.5£). 



Thus, Tph ya {X) in (7.48) is approximated as 



i 

90 



(N 1) ( **b«(X)Tf{X) - l 4c^X)Ml{X)T*{X) 



15 



^b a (X)T a \X) - -fc a (X)Ml(X)T*(X) 



+ 



1 

'12tt 

i 



(N — 1)8 (Ml(X)) Ml(X)T^(X) + 6 (Ml{X)) Ml(X)T^(X) 
(N - 1)6 (-Ml(X)) (-Ml(X)f 2 T^(X) 



+9(-Ml(X)) (-Ml{X)f /2 rW(i) 



(7.60) 
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The moment ( [7.55 ) and (7.56) are approximated, in respective order, as 

dImj phys (X) 



and 



(Q(X)Q(Y)) c = 26 ij 5(X-Y) 



A 



d<p 2 (X) 



<p = o 



"2^ Sa S(X Y) [{N - 1)6 (-M 2 (X)) (-M|(X)) 1/2 T w pQ 
+30(-M 2 (X)) (-M 2 (X)) 1/2 T CT (X) 
- i (TV - 1)0 (-M 2 (X)) M 2 (X) + j (-M 2 (X)) Ml{X) 



<CiWO(y)Cb(tOCi(V0>cc 



A 2 



288?r 



<5 4 (X - CT) <5 4 (I/ - F)<5 4 (X - y) [5 lfc % + 6 kj S H + 8 l3 6, 



lk\ 



x {(N-l)d(-M 2 (X))+96(-M 2 (X)) 



(N-l)9 (-M*(X)) (-MUX))- 1 / 2 %{X) 



-991- 



M 2 (X)) {-Ml{X))-V*%{X)]} 



Im lphys (X) = Jm 7 <° ] (X) + Im~tf h ] (X) vanishes for Ml (X) > and M 2 a {X) > 0, i.e., for 



2 {X) > Mbx(~M^(X), ~M 2 (X)^j 



Then, the condition ( 7.51 ) is met 



2. Low-temperature limit 



We consider locally thermal equilibrium system, which is characterized by ( 7.35 ), and obtain %h ys (X) in the low- 



temperature limit, T(X) << —M 2 (X). The function ( 7.35| ) satisfy the assumption made above after ( 7.43| ), we can 
use the representation Q7.44). From (7.44), one can easily see that the contribution from the segment C( dominates 



over others. Changing the integration variable from z to p through z = (— Mi(X))[p — i], we have 



(0) 



3tp 



(-Ml(X)) 2 / dp 



p i ' 2 {p-2if/ 2 



exp { \P -V)- 1 



-3ni/4 pi (-Ml(X))^ 2 /T(X) 



-M\{X)) 2 [ dp p 3 / 2 exp 
Jo 



2V2 

3^~ £ 

1 e -3../4 e .(^?(^) 1/! /T(i) (-M 2 (X)) 3 / 4 {T{X)f 2 , 

V 27T 



{-M^X)) 1 / 2 
T{X) 



from which, we obtain 



Reif = -^(-Ml(X))^(T(X))^cos ^ 



Im 



7| 



(0) 



V2ir 



{-M 2 AX)f' A {T(X)f 2 



sin 



T(X) 4 

\-M 2 (X)f/ 2 3n 
T(X) 4 



(7.61) 



Incidentally, ( 7.61 ) is directly obtained from ( 7.39|) through an analytic continuation, M 2 (X) (> 0) — * M 2 ^{X) — h 



(M\ < 0). 

The form of l^ ys [X) at low temperature is obtained by substituting ( 7.61 ) into ( 7.48 ). It can be shown that the 
condition (7.51) holds. 
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C. "Determining M%(X) and M 2 a (X) 



In this subsection, following the method 1) in Section VID, we derive the gap equation (the determining equation 
for M^(X) (£ = TTjf)) to the leading one-loop order. Main task is to compute T1r(X;R). Let T>[L\, £2, ■■■] {Cj € 
Ci n t + C rc ) be an one-loop diagram contributing to Xr(X;P), which includes the vertices coming from C\, £2, 
The relevant ones are 



I. V 



II. V 



\(4> ) /4! (tadpole diagram). 
-A(V-<?)0 2 /3!, -A(V-0)0 2 /3! 



Including the contributions coming from C rc in ( |6.10| ), we write Tir{X; R), with obvious notation, in the form 
£ fl (X; R) = (S fl (X ))! + (S K (X; i?)) n + (Z M Z - l)/i 2 

+^(z A z 2 -iv 2 (x)(p.(x) + 3P ff (x))- £ P e POx«PO- 

£=7r, cr 



The leading contribution to ( |7.62| ) takes the form 

V R (X;R)= P{W4 fl (^ fl ) 



£= 7T, cr 



£= 7T, cr 

A 



£<f(X; J2)) n + (Z^Z - l)^ 2 - xe(X) 



V-d^ZxZ 2 ~i)v 2 (x) 

with c? T = 1 and d a — 3. The resummation program completes by imposing the conditions 

ReT:f(X;R 2 = M 2 (X)) = (£ = tt, a). 
This is the gap equation, which serves as determining M 2 {X) (£ = tt, a). 



(7.62) 



(7.63) 



(7.64) 



1. Computation of (En(X)) J 
We start with computing the leading O(X) contribution from the tadpole diagram, (Sr(X))j: 
(Srt=^ | p[2{p,(I)AW(I;P) + P,(I)AW(I;P)} 



+l{(JV-l)A£ ) (X;P) + A£ ) (X;P)} 



(7.65) 



where is the (1,1) component of AW in ( 6.44 ), etc. Note that in ( 6.42 ) yields higher-order contribu- 

tion to Sfl(X), which is pure imaginary and consists of terms that are proportional to (d\if(X))/dX^)(ip(X)\ and 
\(p(X))( d((p(X )\/dX fl ). When one performs two-loop calculation, this contribution should be taken into account. We 
rewrite ( [7.65j ) as 



(Su(x)) I = VAX) [rt*\x)) i + p a (x) [E^(x)) i , 

(xfiX^ = ± [at MX) + /% MX)] (£ = tt, a) 

(K,Ar) = (JV + l, 1), ( a<r , = 1,3)), 



where 



MX) 



d 4 p 



Ai?(X;P). 



(7.66) 
(7.67) 



The O(A) contribution to 2j (X) has already been computed in the above subsections. 
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2. Computation of Re (S_r(X; -R))jj 



Here we compute the leading 0{\ 2 ) contribution to Re (Sjj(X; R)) n - Straightforward manipulation yields 



(E B (X; R)) n = P„{X) R)) ^ + P CT (X) (££>(X; fl)^ 

(4^(X;i?)) n =^(JT)^ <r (Jf;iZ), 

J?)) = ^ 2 {X) [(N - 1)J^(X; R) + %J aa (X- R)\ . 



Here 



d A P r 



(2^) 4 



a£ } (X; P) AiV (X; P - P) - A£ (X; P) A% (X; P - P) 



,(0 



(€), 



,(0/ 



jg\X;R) + jg>(X;R) 



with 



4° 3 (X;P) 



J ? f(X;P) 



(27r) B 



P 2 - M 2 {X) +ir) (R - P) 2 - M 2 (X) + ir\ 



+47T 2 #(-p ) 0( Po - r ) <5(P 2 - X 2 (X)) 5((P - P) 2 - M 2 (X))] 



d 4 P 

(2^)3 



AT(f)(X;|po|;p)5(P 2 -M 2 (X)) 



{R ~ P) 2 - M 2 C (X) - ie( Po ~ r ) V 



Here we have dropped the contributions that have to be dealt with in two- loop order (cf. above after ( 7.6 5| ) ) 



Calculation of 3$ with M^(X), M 2 (X) > is straightforward. The resultant form for may be analytically 
continued to other region than M 2 (X), M 2 (X) > 0: 

e efij 2R Z M~(X) 



+ 



5i 
2P 2 



In 



■M 2 P0-.M 2 p0 + P 2 -S 1 
7W 2 (X)-M2( X ) + il .2 + iSl 



In 



■M 2 (X) - Ml(X)-R 2 +S 1 
MUX) - MUX) - R 2 - S, 

where M^(X) is understood to be M 2 (X) — ie(ro)r] and 

S 1 = V(P 2 - MUX)) 2 - 2MUX)(R 2 + MUX)) + M*(X). 



(7.68) 



Equation (7.68) is valid in the regions, 
Rl) Ml(X)Ml(X) > and 



P 2 < Ml (X) + Ml(X) - 2^MUX)MUX) 



Ml(X) + Ml(X) + 2^Ml(X)MUX) < R 2 . 



R2) M 2 a {X)M 2 UX) < 0. 
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In the region, 



Ml(X) - Ml(X) - \R 2 \ > 0, M 2 a (X), Ml(X) > 0, 



Equation ( |7.68[ ) i s as i t is. The expression in the regions M 2 a {X) - Ml(X) - \R 2 \ < and/or Ml(X), Ml{X) < 
is obtained from (7.65) through simple analytic continuation. 
Continuing (7.68|) to the region 



Ml (X) + Ml(X) - 2^/Ml(X)Ml(X) < R 2 < M 2 a (X) + Ml(X) + 2y/M*(X)Ml(X) 

(Ml(X)Ml(X) > 0), 



we obtain 



(^) 2 7^<x-m- 2 1n ^W , R 2 -Ml(x) + Mj(x) MUX) 

(4.) JiJ(X,R)---ln— + lU Mi(X) 



S2 

R 2 



arctan 



5 2 



arctan 



M 2 a (X) - MUX) + i? 2 

S2 " 

M 2 a {X) - Ml{X) - R 2 



(7.69) 



where 



S 2 = y/2Ml(X)(R 2 + MUX)) - (R 2 - MUX)) 2 - M%{X). 
Equation ( 7.69 ) is valid in the region 

Ml(X) - Ml(X) - \R 2 \ >0. 
If Ml (X) - M\{X) + tR 2 (t = ±) in ( fjxjj ) is negative, the replacement 

1S2 &2 



arctan 



MUX) ~ MUX) + tR 2 



7r — arctan 



MUX) - MUX) - tR 2 



should be made. 

AO) 



Jtf (R) (£ = ^ CT ) is obtained from and ( ggg ) by setting M 2 (A) = -M^(^) = A4| (X). 



As seen in ( |7.64| ), we need iie £^ 5 (AT; .R) at i? 2 = A4|(A). As to J^'(X;R), following [|10[, we compute it 



at 



i? = for simplicity. The difference between J^\X;R = 0) and J^'(X;R) at i? 2 = M'^X) is of higher order. 
Then, we have 



[j^(X;R 2 =M 2 UX))+J^(X;R = Q)_ , 

Re (EW) n ~ ^ ^ 2 (A) [(TV - 1) {j(?(X; i? 2 = X 2 (X)) + J^(X; R = 0)} 
+9 {j-j°> (A; i? 2 = MUX)) + Jg> {X;R = Q)}], 



where 



Jti\X;R 2 
J£ ] (X;R 2 



MUX)) = 
MUX)) = 
MUX)) = 



(4tt) 2 
1 

(I^j 2 
1 

(i^) 2 



•M 2 (A) |M 2 (A)| |M 2 (A)| (T) 
2A4 2 (A) ln |A^ 2 (A)| ' ' '" rK - 

111 O 9 I TT7T 



e 2 ^ 



e 2 ^ 2 



In 



|A4 2 (A)| 



V3 ~ e 2 ^ 2 



(7.70) 



with 
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and 



JCW - 9 (Ml(X)(Ml(X) - 4Ml(X))) 



^Ml(X)(Ml{X)-±Ml{X)) 
2 MUX) 



In 



\Ml(X) - ^Mj{X){Ml{X) - 4MUX))\ 
\Ml(X) + ^Ml(X)(Ml(X)-AMl(X))\ 



In 



\Ml(X) -2Ml(X) - ^Ml(X)(Ml(X) - 4M*(X))\ 
\Ml(X) -2Ml(X) + ^Ml(X)(Ml(X)-4MUX))\ 



(Ml(X)(4Ml(X) - Ml(X))) 



y/Ml(X)(4MUX)-Ml(X)) 
MUX) 



e(Ml(X)) arctan 



y/MUX)(4Ml(X) - Ml(X)) 



-e(Ml(X) - 2Ml(X)) arctan 



\Ml(X)\ 
y/Ml(X)(4Ml(X)-Ml(X)) 



\M*{X)-2Ml(X)\ 
+vr {8 (-MUX)) - 9 (2Ml(X) - Ml(X)) }] , 
9 (Ml(X)(Ml(X) - 4Ml(X))) 



(7.71) 



x ^Ml{X){Ml(X)-AMl(X)) ^ \Ml(X) - ^Ml(X)(Ml(X) - AMl(X))\ 



Ml(X) 

+6 (MUX)(AMl(X) - MUX))) 



\ MUX) + ^MUX)(MUX) 4M*(X))\ 

^/MUX)mUX)-MUX)) 
\Ml(X)\ 



2 arctan 



^MUX){AMl{X)-Ml(X)) 
\MUX)\ 



J^(X;R = Q) 



J^(X;R = 0) 



i ( f\x)-ij p \x) 

MUX)- MUX)' 

dif\x) 

dMj(X) ' 



(7.72) 

(7.73) 
(7.74) 



Here lf\x) is as in ( [Hi]) o r ^Tv^j. E xpressions ( |7J0| ), ( pf77l| ), and ( f7J2| ) do not diverge at Ml(X) = and at 
■M 2 (X) = 0. Equations~^7j|) and ([7.74|) are valid for M\(X) > (£ = 7T, a). Analytic continuation of them to the 
region A4^(X) < may be performed as in subsection B: 



if\x) 



2n 2 



Re 



I dz ^z 2 -M 2 {X){N { ^{X;z)), 



where the integration contour C' is as in (7.44) 



(7.75) 



3. Gap equation 



Let us turn back to the condition (7.64), which may be written as 



Xt = M 2 Ax)-^ 



(X$(X)) I + (4° W) n + (Z,Z l)^ 2 + \d^Z x Z 2 l)^(X) (£ = tt, a). 



(7.76) 



The UV-diverging contributions in (7.76) cancels. In fact, substituting ( 7.19D into ( 7.7C ) , we obtain the UV-divergence 
free equation, 



53 



Here (cf. (7.66)) 



with 



and 



(M|(X)) i + (Mf pf)) n = Mf (X) — n 2 (C = tt, a) 



|A4f(X)| 

-X|(X) In 



(4tt) 2 - e^S 



(7.77) 



(7.78) 



(7.79) 



+Re 



(4n)i\2Ml(X) \Ml(X)\ 

l { f\x)-lj p) (X) 
Ml(X)-MUX) 



In- 



(7- 



(4tt) 2 



(iv-DlmMP-^ 



e 2 /i 2 



e 2 ^ J &M 2 (X) d.M 2 (X) 



(7.81) 



In the above equations, lf ] {X) is as in d7.12| ) for M|(X) > and as in ( |7.75| ) for -M^X) < 0. /C^ and /C^ are as 
in (7.71) and ( 7.72] ), respectively. For the region, 



Ml(X)>0, Ml(X)>4Ml(X) 
or 

M 2 a (X)<0, Ml(X)<AMl(X), 

the first terms on the RHS of ( 7.71 ) and of ( 7.72| ) should be used and the second terms should be used for the region, 

Ml{X) > , Ml (X) < iMl(X) 
or 

M 2 a (X) < , Ml {X) > AMl(X). 



In the LHS of (7.77), we have dropped the higher-order term 

A 1 



48tt 2 e 



(7.82) 



which should be dealt with when one computes next to the leading order. The one- loop diagram 
0X^P^0/2, — A(^> 2 ) 2 /4! , which contributes at next to the leading order, cancels the diverging contribution (fj 



V 



Equation ( 7.77| ) is the gap equation, by which M|(X) (£ = 7r, a) are determined self consistently. M|(X) thus 
determined depend on <p 2 (X). With M|(X) in hand, one can judge if M\{X) in ( |6~2S| ) is positive or negative. 



4- "High-temperature" expansion 



"High-temperature" expansion of lf \x), Eq. (7.75), is obtained as in subsections A and B, 



Relf\X) ~ ^ct(X)T?(X) -±-6(Ml{X)) (M\{X)) 1/2 T c (X). 
From this the "high-temperature" expansion of d2^\x) / dM^(X) (£ = it, <t) in ( 7.81 ) is obtained: 



dif\x) ^ _i_ 

dMj(X) ~ 8vr 



(Mj(x)) {Ml(x))- 1/2 T t {X) 



When M\{X) < 0, dxf \x)/dM\{X) is negligibly small. 



(7.83) 
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5. "Low-temperature" expansion 



In the case of locally thermal equilibrium system at low temperature, we have, for Mi(X) > 0, 



r(/3) 



(X) ~ ^L- (M 2 (X)) 1/4 {T (X)f 2 e- M t^/ T ^ (M 2 (X) > 0). 



(7.84) 



For A4%(X) < 0, computation may be carried out as in subsection B. The result agrees with that obtained from (7.84) 
by an analytic continuation; 



ReX 



03) 



-< w-^(--^iw) 1/4 (W) ^ v T{x) 

From ( |7.84D and ( |7.85| ), the low-temperature expansion of dlf\x)/dM 2 (X) (£ = in ( |7.8l| ) is obtained: 



f(-M £ (X)f 2 r 



(MUX) < 0). 



(7.85) 



Re 



dif\x) i_ 

dM 2 (X) ~ 4^tt 



{M\{X)) (Ml(X)y lfi (T(X)) 1/2 e-MlW/nx) 



-0 (-Ml(X)) sin [ { - M ^f 2 - l) (-MKX))-^ (T(X)f 2 



VIII. APPLICATION — EARLY STAGE OF PHASE TRANSITION 

There are few cases, for which approximate analytic calculations yield reliable results. In this section, for the 
purpose of illustrating how the formalism constructed above works, we apply it to very simple cases, in which the 
system of coupled equations obtained in the previous section can approximately be solved analytically. 



A. Spatially uniform system 

In this subsection, we apply the formalism to a spatially uniform system. We suppose the following situation. At 
the initial time Xq — 0, i) the system is in a "highly super-cooled state" with ip(Xo = 0) ~ and ii) the temperature 
T is still high and the system is in quasi thermal equilibrium, 



aKO( Xo = o ;N) = _±_ 



(8.1) 



We are only interested in the small time interval, during which iii) the high-temperature expansion is applicable, so 
that keeping terms that "explicitly" depend on T is a good approximation (cf. below) and iv) the condensate fields 
ip(Xo) remains small in the sense stated below. As mentioned above, our formalism applies to the system away from 
the critical region, which is the reason to include i) into the list. 

We start with the determining equation for f^(X ;p ) = 9( P o)N^(Xo;po)-9(-po){l+N^(X ; -p )) (£ = vr, a), 
Eq. (|6~53|), 



dN^(X Q ;p ) dN^(X ;p )dM^Xo) 

^PO 1 ~ — — U. 



dX 



dpo 



dX 



(8.2) 



Since N^(X ;p ) in the propagator accompanies S(P 2 — M. 2 (X )), pa in (^2) is po = Jp 2 + A4 2 (X ). Equation 



(3.2) may be solved as 



N^(X ;po) = T{P 2 - M\{X Q )) t 
where T is an arbitrary function. For simplicity of presentation, we take <p(Xq = 0) = d(p{Xo) / dXo 
sequel. In the case of 0{Xq = 0) ^ and/or d<p(Xo) / dXg 



(8.3) 
= in the 



Using (3.1) in (|8.3|), we obtain, 



X = 



X = 

^ 0, no new feature arises, as long as they are small. 
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N^(X ;p a ) = 



1 



e ^/p 2 -M 2 s (X )+Ml(X =0)/T _ l 
1 



= v /p2 +M 2 (0)/T 



(8.4) 



where use has been made of the fact that (Xq;pq) accompanies S(P 2 — M^(Xq)) in the propagator. 

With (^4), let us compute the fundamental quantity T^(Xq), Eq. (7J3). Computation goes in a similar manner as 
the "direct calculation" in Section VIIB: 



^ = 12+2^ 



T 2 



1 



d£ 



dp 



1 



o y Jp l +M*(X ) e y/ p2+M 'W T 



T 2 T 
12 + 2^ 



e £ l T -1 
dp 

o 



^ 2 H>lf(*o)|^-K(0)| 



dp- 



'v^fW ^ 2 - |m|(o; 

where high-T approximation has been used. Noticing that, 

^/p 2 -\M 2 (X Q )\ = ^\M 2 (Xo)\-p 2 for p 2 < |M 2 (X )|, 



(8.5) 



etc., we obtain 



^ 12 



2tt 2 



|M 2 (0)| 



iT 
2^ 2 



^/\M 2 (0)\E[-, 



tt /|M|(Q)|-|M|(X )| 



|M|(0)| 



(8.6) 



where 



iJ(p 2 ) = / dx 



Vl — x 2 \J p 2 — X 2 



dx ■ 



Vx 2 ^! 



r /2 i 

E(%/2,k)= d0Vl-k 2 sm 2 t 
Jo 



Here E is an elliptic function of second kind. Note that |Af 2 (0)| > |.A/^(Xo)|, so that p < 1. From the setup iv) 
above, 

|M 2 (0)| - \Ml(X )\ = Xe^ 2 (X ) « \M 2 (0)\ {e n = 1/6, e CT = 1/2), 



and then (3.6) may be approximated as 

T 2 1 T 



12 4tt 2 
iT 



|M?(0)| 



[|M 2 (0)| - |M|(X )|] 



+ 



4tt 



|A4|(X )| , t 



1|M 2 (0)|-|^ 2 (X )| 



|M|(0)| 



y2 

1~2 



A 



; 4tt 2 



|M|(0)| 



=^ 2 (X ) 



47T 



v^i(*o)i 



Ae 4 ^ 2 (X ) 



|M|(0)| 



(8.7) 
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Let us turn to the determining equation for Mf (X g ) (£ = tt, u), Eq. ( 7.771). As mentioned above, we ignore the 
temperature-independent terms. Using in ( 7.79| ) - ( 7.81 ), we see that ( 7.77 ) yields 



M 2 (X ) ~ -u 2 - 
A 3 



A 
12^ 



(N + 4) 



T 



vmwi 



T 



108^ 2 y/\M$${M*{X ) - Af 2 (X ) + A^ 2 (X )/3' 



M 2 (X ) ~ -c^ 2 - 3 



A 

127T 



where 



(iV + 8) 



A 



T 



<? 2 (*o), 



(8.9) 



^ 72 (iV + 2)T (> 0) - 



In obtaining (8.8) and (S.9), use has been made of the approximation (|Af 2 ( 
justified a posteriori. 

Approximate solution to (3.8) may be written in the form 



|Af 2 (0)|) << |Af 2 (0)|, which will be 



\fz(X fl )^-^-[ A) , V + 1/) 



V\mIW\ 



where 



.9-4 



£ 2 (*o), 



AT 



(8.10) 



VWW) 



For AT < O(V|Af 2 (0)|), g = 8 and, for AT > O(^|Af 2 (0)| ), 5 = 4. 

Let us turn to deriving the equation of motion for <p(X). £l(ip 2 (Xq)) in (7.24) reads 



n(0 2 (x Q ))~~uj 2 + ^ 2 (Xo) 

6 



AT 
"12^2 

Up to 0(</? 2 ) term, this becomes 
where 



(N-l)y/\Mm\H 



f \Ml(X )\ 
V \M 2 (0)\ 



■3^\Mm\H 



( \Ml(X )\ 
V |M 2 (0)| 



n(^ 2 (X )) ~ -uj 2 + ^0 2 (Xo), 
b 



A' = A 



247r 2 v ' u 



.11) 



It should be emphasized that A' < A. Above equations ar e val id for AT < 0(uj) (cf. the setup i) at the beginning of 
this subsection). The equation of motion for 0{Xq), Eq. ( 7.54 ), reads 



where <p(Xo) = d 2 tp(X ) / dX^ and 



A' 



0{X O ) + C(X ) + He 



dV 



difiiXo] 



A' 



V ~ -NT* - 2 (X a ) + - {0 2 {X O )Y HCe- 0(X O )) 



4! 



.12) 
.13) 

.14) 



57 



is related to Jphys through 



V 



-Rej p hy 



^Vo) 



with (p{X(j) = dip(Xo)/dXo. In the following we simply call V the potential. From the setup i) at the beginning of 
this section, u 2 ~ — M 2 (0) ~ — M%(0) (> 0) is quite large. The setup iv) means that we are interested in the time 
interval during which oj 2 » \'0 2 {Xq). [Recalling that X of <p(X ) is macroscopic time coordinate, uo should be 
oj < 1/io, where Lq is the size of the time-direction of the spacetime cell (cf. Section I).] Above observation allows 
us to solve ( g. 12 ) iteratively. Ignoring the term \'0 2 (Xo)<p(Xo)/6, we obtain the zeroth-order solution, <p^(Xo), of 



(X ) = ^ [cosh(iuX Q ) - 1] g+ - I dX' ({X' ) sinh(a;(Xo - X' Q )) 



H 



Xn 



.15) 



Substituting ( S.15 ) back into the term AV 2 (^o)^(^o)/6 in ( 8.12| ), we obtain the first-order solution, 



<p(X Q ) = —r [cosh(wXo) - 1] e 



dX' a 



ax' )-^{^\x' Q )Y^{x' Q ) 



sinh(^(X - X' )), 



(8.16) 



where £ (0) (^o) is as in ( |B.15| ). More accurate solution may be obtained by further iterations. 

For the purpose of computing the ^-averaged correlation functions, we compute the partial derivatives of 
ImjphysiXo) with respect to ip 2 (X ) using (7.24) and (7.25) with (jB7?|): 



d Im"f phys (X ) 



d0 2 (X o ) 
d 2 Im^physjXp) 



0=0 487T 

^^(TV + 8)^. 

ip = 11527T LU 



(8.17) 
(8.18) 



d& 2 (X )) 2 

Noticing that we are concerned about the spatially uniform system, (7.55) and ( 7.56| ) become, in respective order, 



and 



UX a )Q(Y )) c ~ -±-{N + 2)^ 5 l3 S(X a - Y Q ) 



\ 2 iu I o T 

<C 4 (x )o(r )a(c/o)0(K)))cc ^ - 2887r ya ~ Wn + M« + Wik] 

xS(X - Y Q ) 5(X - U Q ) 5(V - Y ), 



(8.19) 



(8.20) 



where V is the (large) volume of the system. 

We are now in a position to co mpute the ^-averaged "connected" correlation functions ( 7.58j ). From the zeroth-order 
solution ( B.15 ), we obtain using ( S.19 ) and ( 8.20 ), 



(v? (0) (^o)> C 

(^ 0) (x )^ 0) (ro))cc 



h 



e [cosh(wXo) — 1] 



.21) 



— y 



pae 



Lo(pX -aY a ) 



dX' 



A T Si 



dY'e-^x'o-^^X'MY^ 



- — . -£(N + 2) [cosh(u;X ) sinh(u,Yo) - u Y cosh(c i ;(X -Y ))] 

487T w z V 



.22) 



and 
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(^ \x )^>(y )^"(u q )^"(Vo)) Cc 



JO) 



(0), 



[SijSki + 8 ik Sji + SuSjk] 



N + 8 ( A V T 

"hT vis/ ^v" 3 

x [cosh(w(X + Y Q + U Q - V )) sinh(2wVo) 

- Asmh(LuV ) {cosh(w(X + Y - V )) cosh(w([/ - V )) 

+ cosh(w(X — cosh(wC/o)} 

+ 2wV {2 cosh(w(X - Vq)) cosh(w(r - [/„)) + cosh(u;(X -Y -U Q + V ))}] 



In ( 3.22j ), Xq > lo has been assumed and in ( 8.23| ), Xq, Yq, Uq > Vq has been assumed. Let us study ( |8.21 ) 
in turn. 



.23) 

H) 



1- & {0) (x )) c 



Equation ( |S.21 ) shows that, as Xq increases, the £-averaged tp(°\Xo) blows up exponentially, 

(^ \x Q )) c ~$e" x <> (uX »l), 



.24) 



where $ = He/(2uj 2 ). We should recall here that, for too large Xq, 0(Xq) is large and our app roxim ation does not 
apply. For obtaining an improved expression for (0(Xo))^, we should use the improved solution (8.16), which yields, 
for large e ujXo » 1, 



<«>)>< 



uiX 



A' 



48w 2 



a ,„ ,., a ir 



1927T 



(iV + 2) — — 

W 2 K 



.25) 



Comparing ( jjTgg ) with ( gTgg ), we see that Q3.25j ) is reliable in the region 



1927T 



W 2 ^ 



which generally forces A to small values. When compared to the first term on the RHS of ( B.25| ), the term with $ 2 
in the second term prevents |(<p(Xo))| to become large, while the other term expedites it. The former can easily be 
understood from ( B.12| ) , since the term \'0^ 2 (X o )0 w (X o )/6 in (§3|) acts as decreasing ui. The latter comes from 



~ {<P (a) (Xo))\ ( ?\x a ) -> ^((^(X ))\^\X ) + (X )^°\x ) kip f\Xo) 



(8.26) 



in (3.12). As seen from (8.22), the ^-averaged quantities in (8.26) are negative so that, in (8.12), ( 8.26Q acts as 
increasing |<p(Xo)| as Xq increases. 



2. Two-point (-correlation function 



We now turn to ( 8.22| ). One can easily see that flS.22 ) is negative definite. Thus, Im^ p hy S {Xo) causes a negative 
correlation between the same component of ip, <pi {Xq) and tpj = i(Yo). If the term ((X ) in ( |3.12| ), which comes from 



I^njphys(Xo), is absent, (8.15) or (3. If) with ( = uniquely determines ip(Xq). Imj p hys works as stochastically 
upsetting <p(Xo) through the random forces <f. Since ((p\ ' (Xo)ipl°\Yo))^ (with no summation over i) in ( |8.1S| ) is 
negative, (ipi(Xo)ipj—i(Yo))( t tends to decrease when compared to the case where no random force acts, C(-Xo) = 0. 



As Xq increases, the £-averaged two-point correlation function (8.22) blows up exponentially. For X = Y and 



» 1, we have 



{ipf\x Q )yf\x Q )) 



1927T 



(N + 2)^ S 4e 2 " x °. 

LO z V 



.27) 
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Using the improved solution (8.16), we may compute the im prov ed (- averaged two-point correlation function. Up to 
the first order with respect to the difference between flS.16 ) and flS.15 ), we obtain 



(^(Xo)^(X )) Cc ~ Eq. (8.27) +C^(X ) + C {2) (X ), 



where 



,(i) 



,(2) 



AA' 



1 T „4. 



4608^ 1 + Vu 



da® 2 + 2<f> t $ 



%A 2 A' 



N ■ 



96 



1 T 



»3* ' "^1^3 
2 



1 - 2tt 



iV + 8 



(N + 2) 2 w 2 TU 



(8.28) 

(8.29) 
(8.30) 



Let us write ipi(Xo), Eq. (8.16), in the form 



i p i (X )=Y l cpf(X ), 



1=0 



where ipf ] (X ) includes I C's. c}j 3 (X ), Eq. ( gTgg ), comes from (^f 1 (X )^ 1] (X )) c and C^ 2) (X ) comes from 
(v[ (Xo)v^ (-^o))c and (y|"(X )^ (X ))^. Q- (Xo) consists of two terms. The term being proportional to 1/V 2 
comes from (the square of) (8.19) while the term be ing p rop ortio nal to 1/U 3 comes from (8.20). 



Here we summarize the characteristic features of ( 8.28 ) - ( |B.3C| ) 



• Comparing ( |8.28j ) - flS~3C| ) with ( |8.27[ ) we see that ( |8.28j ) with flS.29| ) and ( |S.30D is reliable in the region, 

2 LJ 4 V LJ 6 V 2 ' 



e 2uX < q Min 



A'$ 2 ' AA'T' AA' 



• In contrast to flS.27 ) and C^\Xq), C^\Xq) is not diagonal in TV-dimensional vector space, which means that 
the correlation between different component of ^(^o) exists. cjp(Xo) (with no summation over i) is positive, 
while C^ i} (X ) can either be positive or negative depending on 



The 1/V 2 term in C\, (Xo) is negative and then enhances (8.27), while the 1/U 3 term is positive. 



3. Four-point (^-correlation function 
Setting X Q = Y Q = U = U in JO^) , we obtain 

<^ 0) (x )^ 0) (y )^ 0) (c/ )^ 0) (u )) Cc 

[SijSki + S ik Sji + SaSjk] H(ujX ), (8.31) 



3 

N + 8 / X\ 2 T 



4vr V 48 / uj 6 V 3 
where 

H{ujX q ) = sinh(4wX ) - 8sinh(2wX ) + 12wXi 



One can easily show that TL{loXq) (X > 0) is positive definite, so that the four-point correlation (8.31) is of negative. 
For uiXq << 1 and uiXq >> 1, we have 



n(ujx ) 



f(u,X f (u,X « 1) 
\e iuX ° (ljX » 1). 
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4- Features of the distribution function (8.4) 
For the purpose of comparison, instead of ( |S.4D , we take the quasi equilibrium distribution function, 

N^(X ; Po ) 



rcf 



>y /p 2 -\M 2 ( (X )\/l 



Same approximations as above yield, 



12 4tt 



yJ\M}(X )\, 

0{X O ) + C(X ) + He, 



dV 



d0(X o ] 



A .(iv+2: CJ(i ' r 

A 2 N + 8 T 



1447T lu 

x5{X - Y ) S(X - U ) S(V - Y Q ). 



Let us summarize distinctive features of our results when compared to (8.32) - (3.38). 

• Change in A, A — > A' [Eq. (8.11)], does not occur in (HH) - flOil ). We emphasize again that A' < A. 



• We rewrite the potential V, Eq. (8.14), in the form 



Here 



w 2 + Aw 2 = oj 2 + 2(N - 
A 



A \ T 



24tt / w 



¥? 2 (*o) 



with 



AT 



V-^(^ + 2)(T + AT) 2 



A iV + 8^ 2 (X ) 



8tt 2 N + 2 w 



Equation (B.39) is to be compared with (8.36). Comparing ur + Aw with 



J = -f-±{N + 2)T\ 



(8.32) 

(8.33) 
(8.34) 
(8.35) 

(8.36) 
(8.37) 

(8.38) 



Equa tion ( ft.38j ) is two times ( ft.20| ), so that, besides the change A -> A', ( |8.23| ), l/V^part of Q ^ in Q8.30| ), and 
(3.31) arc one-half of the corresponding formulae with the reference distribution ( |3.32] ). 



(8.39) 



(8.40) 



(8.41) 



we see that, as <f(Xo) is "rolling down" the potential hill, the temperature "decreases," and thus the system 
"cools." 
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• We write the distribution function (8.4) in the form, 



N^\X ;po) = 



1 



eV /p2_|JW2 (Xo) |_ esA v?2( X() ) /T _ x 
1 



^p>-\M\{X Q )\/T eff (X ) _ 1 



(e w = 1/6, e ff = 1/2), 



.42) 



where 



T 



p 2 -\M 2 (X o )\-e i X0 2 (X o ) 



(8.43) 



Let us co mpare this with the reference distribution (8.32). We hrst note that the parts of the stable (oscillator) 
modes in (8.32), 



\M}(X )\ <p 2 < \M 2 (X )\+e s \tp 2 (X ) 



(8.44) 



are the unstable modes in ( 8.42j) and thus the "number" of unstable modes in the latter are larger than that in 
the former. For the modes, pVt|(Xo)| + e^Xifi 2 (X ) > p 2 , ( 8.43 ) shows that the effective temperature increases 
when compared to the reference distribution ( |3.32j ). This temperature increase is interpreted as follows. As 
seen above, while <p(Xo) becomes large, the potential V becomes small and the energy of the condensate fields 
0{Xq) is released and converted to th e kinetic energy of the quantum fields. On the other hand, for the 
unstable modes, p 2 < \M^(Xq)\, ( 8.43 ) shows that the effective temperature decreases when compared to the 



reference distribution fl8.32Q . These unstable modes yield negative contribution to ReT^ in (3.6). All the effects 



mentioned above compete and ReXi^ turns out to be smaller than the reference-distribution counterpart. 



B. Nonuniform system 

In this subsection, we briefly analyze the spatially nonuniform system under the similar situation as in the last 



subsection. We assume that the initial distribution function is given by 

1 



.1 



N®(X =0,X;P) = 



Apo\/t _i 



(8.45) 



The determining equation for f®(X;P) = 6(p )N& (X; P) - 6(-p )(l + N®(X; -P)), Eq. §1% , may be solved 
as in (|1|): 

N®(X; P) ~ F(X- P 2 - M\{X), p), (8.46) 

where X and p are as in ( 5.43| ) with T^„ = 0. Substitution of the initial data ( 8.45| ) into ( 8.46| ) shows that T is 
independent of X and p and N^(X; P) becomes (|S~4|): 

N^(X;P) = 



>v /p 2 +M|(X =0,X)/2 



At the initial time Xq = 0, \ip\ and \d0/dX^\ are small, so that (cf. (8.1C) and ( |8.9| )) 

M|(A = 0,X)~-^ 2 ^ = 7r,ff), 
which is independent of X. In place of ( 8.1 9| ) and ( 8.20| ), we have, in respective order, 

(CipOCi(*0>C - -^(N + 2) Tu8 tj 8\X - Y), 



.47) 



UX)Q(Y)UUMV)) Cc 



X 2 T 
~^^( N + 8 )- IMfe* + Safin + SaSjk] 

Z007T LU 

xS 4 (X - Y)S 4 (X - U)5 4 {V - Y). 
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In this subsection, we are only concerned with a zeroth-order solution for (p{X) (cf. (8.15)) 



$W(X) = $W e {X)+ I j^—r I dX' o aX' ;p)sinh(cj\X -X^) 7 (8.49) 



(2tt) 3 uj 1 



where 



(8.50) 
(8.51) 



J = ^Juj 2 - p 2 , 

C(x ;p) = / d 3 xe-^c(xix). 

For, lo 2 — p 2 < 0, to' should be understood to be u>' = iy/p 2 — uj 2 . In (8.49), <fipure(X) is the solution to the equation 

(al-c 2 )<4 u ) re (X) = i/e, 

under a given initial data. 



1. Two-point (^-correlation function 
Using ( p7| ), one can readily compute the two-point ^-correlation function 

(<pl 0) (x)<p?> (y)) Cc =-_ a„,.y-: 2u-t 



A ~ _ /" d 3 p jp.(x-Y) J_ 



96tt 



(2tt) 3 



2 2 
UT — £T 



sinh(y/^ 2 - p 2 (X + y )) - sinh^ 2 - p 2 (X - Yq) 
\] io 2 — p 2 



-2F cosh(V^ 2 -p 2 (^o - io)) 

^ f shlh^ 2 - p 2 (X + Yq)) - Sinh(y/^ 2 - p 2 (X - Yq)) 

[ ^u 2 -p 2 

ocosh(^ 2 -p 2 (X -y ))} 



-2Y r 



+9{p 2 - u 2 ) 



sin 



(VF r3 ^(*o + Yb)) - sin(v^^(X - Yo)) 



\Jp 2 — uj 2 



-2Y cos(^p 2 -uj 2 (X - Y ))} 
where Xq > Yq has been assumed. Manipulation of ( 8. 52] ) in Appendix E yields 

(<p<f\x)<pf\Y)) (c = -j^(N + 2)6 l3 



192tt 3V " ' '"' J |X-Y| 
x [6(X + Y - |X - Y|) 0(|X - Y| - (X - Y j) 

x {^(X +Yo)cos(w|X-Y|)-sui(a;|X-Y|) 
l 



sin(u;|X - Y|0 cosh Ul - ?l}(X + Y )) 



-sin(u,|X-Y|)]} 
+0(X O - Y - |X - Y|) 
/ r 1 ^sm(^pC-YJ0 



J2 ^cosh (Vl -£MX + rY ) 



-2wYo / d£ 
'0 



1 J ^sin(a;|X-Y|0 . 



1-e 2 



sinh 



(v / W 2 co(Xo - Yo)) 



(8.52) 



(8.53) 
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First of all, we note that (ip^ (X)^'' (Y))^ c vanishes for |X — Y| > Xq + Yo, which means that the correlation does 
not spread with super- light velocity. Several limiting cases are analyzed in Appendix E. 
1) X Q - Y < |X - Y| < X + Y ; uj\X - Y| < lu(X + Y ) « 1: 

(<Pi{X)<Pi(Y)) Cc * (N + 2)<%- - Cl 



jo) 



384tt 2 
*o+Y 



IX -Y| 



IX-YI 1 



u 2 (X + Y ) 2 



(8.54) 



This shows that the correlation is negative. For fixed |X| and |Y|, the correlation is maximum at the smallest |X — Y|. 
|X - Y| = Max(X - Yo, ||X| - |Y||) and is minimum at the largest |X - Y|, |X - Y| = Mm(X Q + Y , |X| + |Y|). 
We are dealing with the zeroth-order solution to the evolution equation (the counterpart to (8.12)), i.e., the nonlinear 
term, —X'ip 2 (p/6, in the evolution equation has been ignored. The condition that the ignored term is small restrict 
the domain, where (8.54) is valid, to 



AA ,X + Yo- |X-Y| < Q 



|X — Y 

2) 1 << uj{X + Yo) and w|X - Y| << ^/uj{X + Y ): 

/2?r 



(^(X)^-(F)) 



768tt 3 



X(N + 2)<%wT 



1 



(u(X + Y ))3/2 

1 



exp [lj < X + Y 



|X- Y| 2 

2(X + Y ) 



x 1 + 



M x o + Y ) 

+V2^9(X Q + Yo - |X - Y|)0(|X - Y| - (Xq - Y )) 
Xo + Yo 



IX- Y 



cosHX - Y| 



(8.55) 



which is valid for 



XX > e ^x 0+ Y ) < (Y. {lu{Xq + y o)} 3/2 



a n' 1 » + 7 « .^( w ^ 
and AA — — < O {-) 



Equation ( |3.55 ) indicates that, for large Xo +Yo with o>|X — Y| << y/uj(Xo + Yo), the (negative) correlation spreads 
over the distance, 



IX-YI 



2(X + Y ) 



and asymptotically spreads over whole space region. It should be emphasized, however, that, for too large Xq + Yo, 
our whole approximations in this section cease to hold. 
3) |X- Y| < X - Yo; lo(X + Y ) « 1: 



(Y)) C c 



A 

" 384tt 2 



(N + 2)S tj ojT(wY f 



2. Four-point (^-correlation function 



Using (3.49) and (8.48), we obtain 



A 2 T 
~ ~ 2887T^ + ^ uj ^ <li2< ^ 3 * 4 + hiia&iiu + ^114^213] 



Xn 



dX'o j dX! Yl 



Pi 



»Pi'( x 3-x') 



(2tt) 3 



IjS 1 - p 2 



sinh 



(^/u 2 - p 2 (X j0 - X>) 



.56) 



04 



where Ai ,A 2 o,A 3 o > X 40 has been assumed. For equal times, X w = X 2 q = X 30 = X 40 (= X ), computation of 
( 3.56 ) in the limit, ujXq >> 1 and |Xj — Xy| << Xq — 1,..,4) is relatively simple. Similar calculation as in 
Appendix E yields 



(<pf? (X , X x )^f (X , X 2 )^ 0) (X 0) X 3 )^ ( 4 0) (A , X 4 )) c 
A 2 

' (N + 8) [6 



288tt 



(8toI ) 9/2 \ 8A 



which is valid for 



It should be emphasized that, if we use the reference distribution (|8.32|), we obtain two times (8.57) 



IX. CONCLUDING REMARKS 



In this paper, two related subjects are dealt with: 

A) Deduction from first principles of the perturbative framework for dealing with out-of-equilibrium relativistic 

complex-scalar-field systems. 

B) Deduction from first principles of the perturbative-loop-expansion framework for studying an O(N) linear a 

model. 

We have assumed the existence of two different spacetime scales, the microscopic and macroscopic. The first small 
scale characterizes the microscopic correlations and the second large scale is inherent in the relaxation phenomena. 
Besides this setup, in principle, we have not made use of any further approximations. Let us summarize what has 
been added and what has been clarified in this paper. Some additional comments are also given. 



3. On the issue A) 

(1) The assumption that is necessary for the propagator to take the standard form has been singled out. In most 



available work, among 2n-point initial correlation functions W 2rl (n > 1), Eq. (3.13), ( n > 2) are dropped 
and only W2 is kept into the propagator. This is so done as an approximation or by assuming some specific 
form for the density matrix. In this paper, we have seen that, to what extent, W2n (n > 2) may be discarded. 
Whenever necessary, one can incorporate them into the perturbative framework (cf. Q). 

(2) We have emphasized that the perturbative framework is effective only when the inverse size of a spacetime cell 

is much smaller than the infrared scale in the microscopic sector of the theory. In the present complex-scalar- 
field theory, this condition is (L^)^ 1 < \f\V (/i = 0,1,2,3) with A the coupling constant and V the typical 
parameter(s) (of mass dimension) of the system. 

(3) We have proposed mutually equivalent two perturbative schemes. Both of them lead to a generalized relativistic 

kinetic or Boltzmann equation. Derivation of it has proceeded in a physically transparent manner. In the 
physical-/ scheme, the equation is derived through the process of redefinition of the number-density or relativistic 
Wigner function, from the initial number-density function to the physically sensible one. Traditional derivation of 
the equation proceeds (|Jl^l in a rather abstract way, using the Dyson equation. (For the case of nonrelativistic 
many-body theory, see, e.g., p9[|.) As seen in Section V, the Dyson equation is the equation that serves as 
determining the full propagator in terms of the self-energy part. In other word, the Dyson equation simply 
serves as an efficient way of resumming the self-energy part. Then, by referring to the full propagator thus 
obtained, one introduces p| fl3|]39| ] a redefined distribution function. It is this function that subjects to the 
relativistic Boltzmann equation. 
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(4) In massless theories, divergence due to infrared and/or mass (or collinear) singularities appear in some amplitudes. 
This occurs as a result of an interplay of b are massless propagators. By performing resummations of the self- 
energy part for such propagators (cf. ( |5.11 )), one can get rid of this divergence disaster (cf., e.g., p||io||. For a 
different calculational scheme, in which no divergence emerges, see Pq|. 



4- On the issue B) 



(1) Through introducing effective masses, which depend on the macroscopic spacetime coordinates, we have derived 

self-consistent gap equation for them. Together with the generalized relativistic Boltzmann equation (cf. item 
(3) above) and the equation of motion for the condensate fields, this constitutes the system of coupled equations, 
which describes how the phase transition proceeds. It should be emphasized, however, that, as for other work 
of this sort, the present scheme applies only for the systems away from the critical region. 

(2) The negative curvature region of the "potential" has been dealt with by introducing random-force fields. To 

leading one-loop order, this leads to negative correlation between the condensate fields. 

(3) As in the issue A) above, within the gradient approximation, the propagator consists of two parts, A = A.^ + 

AW, where A.^ is the dominant part and the "correction part" A^ 1 - 1 contains a derivative with respect to the 
macroscopic spacetime coordinates. In a (7?, cr)-space, A^ ^ is diagonal, while A^ 1 ) is not. This means that, 
in the level of A^, there occurs mixing between 7? and a sectors. Then, when one takes the part A^ 1 ) into 
account, rediagonalization of the propagator in the (tt, <r)-space is necessary, which leads to "physical" 7? and a. 

Finally, we mention two related works QJl^,^]. In the theories developed there, a "generalized mass (function)" 
is introduced, which is 2 x 2 matrix in the "complex-time plane" and matrix elements are complex functions of 
spacetime. Due to the last fact, the theories provide for a treatment of dissipative effects. The "generalized mass 
(function)" is determined self consistently through renormalization [^,^) or through some optimization procedure 
pOfl . Through this process, the Boltzmann equation emerges. As applied to the complex-scalar fiel d the or y, the part 
of the hat-Lagrangian density that include the "generalized mass (function)" corresponds to (cf. (4.39) - (4.41) and 
(p6|)) 



(with C c as in (4.41)) in this paper. 



^2(-y~ 1 ^(x)M 2 (x)^(x)-C c 



3 = 1 
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APPENDIX A: MASS-DERIVATIVE FORMULA 



Here an argument for necessity of adopting the \po | prescription is given. We add C r 
density fl3.3|). The effect of C rn on the propagator A (A; P) should be 



A(X;P)^ A(X;P) 



m — * m + 5m 



to the Lagrangian 



(Al) 



Let us evaluate the 0(Sm 2 ) correction, 5A, to A due to the perturbation £ r , 



5A(x,y)=Sm 2 / d 4 u A(x, u) T3 A(u, y) 



dm Id u 



XT 3 



A(X,x-u) 



u-y 



d x A{X,x~u) 



A II — T 

A(X, u-y) + -— ■ d x A{X, u-y) 



GG 



where X = (x + y)/2. Taking the Fourier transform with respect to x — y and using A + f^A + = and (4.33), we 
obtain 



6A(X; P) ~ Sm 2 A(X; P) f 3 A(X; P). 



(A2) 



It is to be noted tha t (A2) is valid to the gradient approximation, i.e., the term with first derivative (with respect to 
X) is absent in ( A2). 

Substituting (4.28), we finally obtain 



SA(X;P) ~ Sm 2 



A 2 
A 2 - A 2 



f(A 2 R -A 2 A )A, 



(A3) 



If we adopt the \p \ prescription, / = f(X;P) is independent of m 2 (cf. (4.30)), so that SA in (A3) may be written 
as 



5A(X;P) = Sm 



2 dA(X;P) 



dm 2 



(A4) 



which is in accord with (Al). On the other hand, if we take original / (Eq. ( 4.30| )), / = 0(pa)N + (X; E p , p) — 6(—p )[l + 
N + (X; E p , — p)], it depends on m 2 through E p and 6 A in (|A3|) cannot be written in the "mass-derivative form" (|A~|). 
We have confirmed the consistency between the mass-derivative formula and the \po\ prescription at least up to the 
terms with second-order derivative with respect to X . 



APPENDIX B: PROPERTIES OF THE SELF-ENERGY PART 

B.l. Proof of En + E 12 + £ 2 i + E 22 = 

Let T*i(x, y; Z\, z n ) be a contribution from a loop diagram to the self-energy part in configuration space, where x 
and y are the space-time coordinates of external vertices and zi, ...,z n are those of internal vertices. Here an external 
vertex means a vertex to which an external leg is attached. The times Xo,yo, z±o, z n Q are on the closed time path 
Ci©C 2 . 

Let us first show that, when one of the internal vertices, say Zio, is the largest time, 



z i0 > x , y , other z 's, 



(Bl) 



vanishes Q. This can be proved as follows. From fl3.4| ) and ( |3.5p , we see that a vertex factor with z^ £ C\ 
possesses an opposite sign relative to the corresponding vertex factor with z^o € Ci : 



(iV 4 )i 



-(iVih, (iV 2 )i = -(iV 2 ) 2 , {iVl) 1 = -{iVl) 2 . 



(B2) 



On the other hand, as can be easily be seen from fl2.3|) , all the propagators in are the same for z^ £ C\ and 
Let (S;)«j be T*i(x, y; z%, z n ) with xq £ Ci and yo £ Cj. Similar observation as above shows that 

(Ei)u + (S;) 2i = (i = l,2) for x > y ,Zi,...,z n , 
(S;)ii + (Ei)i2 = (i = 1, 2) for y > x ,Zi,...,z n . 

Then, we obtain 



X] C^i{x,y;z 1 ,...,z n )) lJ = 0. 



(B3) 



There are contributions to the self-energy part, E;(x; Zi, z n ) which include only one external vertex. In such cases, 
(E;)i2 = (S;)2i = 0. Similar argument as above shows that 



^(E ; (x; zx, —, z n ))i 



(B4) 
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There is an exception to this argument, i.e., a contribution from the lowest-order tadpole diagram, (S;(x))y . Explicit 
computation shows that Y^i=i{^l{ x ))ii = f° r the complex-scalar field theory analyzed in Sections III-V and for the 
O(N) linear a model in the the symmetric phase (cf. Section VIIA). This is not the case for the broken-symmetric 
phase (cf. Section VIIB). As a matter of fact, in the region, where the curvature of the "potential" is negative, 

imaginary part develops in En and £22: im£n = ImY> 22 . Then J2i=i(^i( x ))a 7^ but 12i=i R e (^i( x ))u — 0- Thus, 
in studying this region, one should take this f act into account. 

Inclusion of the vertices iV c {x, y) (Eq. ( 4.42 )) that comes from the counter Lagrangian does not invalidate the above 
argument. The vertex iV c (x,y) yields a contribution to the self-energy part Ti c (x,y) = —V c (x,y) (see (5.1)). From 
(fL42|) with J3.16D , it is obvious that 



( E c(x,y))y = 0. 



iV c also app ears as an internal vert ex of S;'s 
above after (Bl) holds, so that 
This completes the proof. 



) and (B4) hold 



Using the form (4.42), we can readily see that the similar argument as 



B.2. Proof of T, A (X;P) = [T. R (X; P)]" 

From the definition of G's, Eq. ([2~g|), it can easily be shown that 

iG n (x,y) = [iG 22 (y,x)]* , 
«G 12 (2i) (x, y) = [iG 12 (2i)(v,x)]* . 



(B5) 



This shows that, by taking a complex conjugate of iGji(x,y), the index j [I] at the end point x [y] changes; j = 1 
(respect. 2) — > 2 (respect. 1) [1 = 1 (respect. 2) — * 2 (respect. 1)]. Let us define 



^i(x,y) 



J ^n^ 4 ^'^ ^i(x,y- 1 z 1 ,...,z n ), 



where E;(x, y; Z\, z n ) is as in the previous subsection. From (3^4) and (3.5), we obtain 



i(v i ) 2 = [i(v 4 ) 1 y, i(y 2 ) 3 = [i(y 2 ) 1 ]*, i{vi) 2 = \i{vi) 



ij • 



(B6) 



Thus, again by taking a complex conjugate of HV^jj, etc., the index j changes; j = 1 (respect 2) — » 2 (respect 1). For 
the O(N) linear a model, similar relations hold. As to the vertex iV c (x, y), Eq. (4.42), we have iV c (x, y) = {iV c (y, x)}* , 
so that 



i{V c {x,y)) 11 = {i{V c {y,x)) 22 }*, 
i(V c (x,y)) 12(21) = {i(V c (y,a;))i2(2i)}*. 



(B7) 



Incidentally, this shows that (T, c (x,y))ij (i, j = 1,2) in (5.1) satisfies (B7) with (V c )ij — > (£c)r/- 
Using (B5) - (B7), we can show that 

- i(Ei(x, y)) n = H(^(y, x)) 22 \* , 
-i(^i(x,y)) 12(21) = [-i(Zi(y,x)) 12{21) ]* . 

Fourier transformation on x — y yields 

(S,(A; P)) n = -(£,(*; P))* 22 , E n (A; P) = -E* 22 (X; P), 

(E;(A;P))i2 , (£;(X;P)) 2 i : pure imaginary. 



(B8) 



(B9) 



Using this in (Oh, we obtain XU(A; P) = [T, R {X; P)]" 
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APPENDIX C: USEFUL FORMULA 



We deal with 



d 4 (x-y)e lP < x - y) j d 4 ud 4 vA(x,u)B{u,v)C(v,y). (CI) 

We wish to express this quantity using the macroscopic coordinate X — (x + y)/2 up to the first order with respect 
to the X dependence. To this end, we expand A as 

A(x, u) = A((x + u)/2; x-u) 

~ A(X; x-u) + — ^ • d x A(X; x - u). 



B and C may be expanded similarly. Substituting these into (|Cl| ) an d Fourier transforming with respect to the 
relative coordinates, we obtain 



dA dB „ A dB dC dA dBC dAB dC 
-C - A 



dP» dX u 



dX^ dP>* dX^ dP» dP» dX^ 



Eq. (CI) ~ ABC+- 

where A = A{X\ P), etc. 

The second form we deal with is 



d i {x-y)e lP - {x - y) j d i ud i vf{x)A(x,u)B{u 7 v)C(v,y)g{y) 
where f(x) and g(y) depend weakly on x. In a similar manner as above, we obtain 

fg 

where / = f(X) and g = g(X). 



Eq. (C3)~fgABC+- 



dA dB „ , dB dC dA dBC dAB dC 

fj A . | 

dP^ dX u 



dXf, dP* 1 dXf, dP>* dP>* dX^ 



dABC 



(C2) 



(C3) 



(C4) 



APPENDIX D: RENORMALIZ ABILITY OF OF THE THEORY BASED ON (5.46) 

Let us call the perturbation theory based on ( |3.3D [( |5.46| )] the theory I [II]. We first observe that, in theory II, 
if we include the term x^x)^ 4> = (M 2 (x) — m 2 )^ ' <j> in C rc into Co, theory II reduces to theory I, which is UV- 
renormalizable. As stated in Section VE, the difference between theory I and theory II lies in the fact that, in theory 
II, x( x )^ $ i s regarded as a perturbative term. 

Let 

n 

T(m 2 ,X) = J2^ l) (m 2 ,X) (Dl) 

i=i„ 



be a renormalized and then finite Green function evaluated in theory I. In (Df ), X is the center-of-mass coordinates of 
the "reaction" described by T and irrelevant arguments have been dropped. If we are making a traditional perturbative 
expansion, stands for the 0(X L ) contribution and, if we are making a loop or S expansion, stands for the 
/-loop contribution. The X dependence of T^'(m 2 , X) has come from various sources. Among those, the relevant one 
to the following argument is the distribution function (in theory I), which we write f\{X;P). In fact, T^ l >{rn 2 ,X) is 
written as a sum of integrals over momenta P's, whose integrands contain fi(X; P), dfi(X; P)/dX ll , etc. In order to 



visualize this dependence, we write (Dl) as 



T(m 2 ,{fi(X)}) = J2^ l Hm 2 Afi(X)}) (D2) 



l=l 

Substituting 
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into rW(m,{/i(X)}), we obtaii 



m 2 = M 2 (X) - X (X), 



t( 1 \m 2 (x)~x(x), {MX)}). 



(D3) 



As stated in Section VE, in theory II, x( x ) is one-order higher than M 2 (x). More precisely, in traditional perturbative 
expansion, x( x ) — O(X) and, in loop expansion, x( x ) = 0(6). Then, we write (|D3j) as 



Eq. (D3) = r«(M 2 (X), {fj(X)}) + £(-*(*)) 



i'=i 



, d l 'T( l \M 2 (X),{MX)}) 

d(AP(x)y 



(D4) 



Note that (-x) l 'd 1 ' /d(M 2 ) 1 ' is of (I + l')th order. 

The distribution functions in theory I and in theory II are different. The UV-divergence issue is foreign to them. 
We can rewrite (M 2 (X), {MX)}) in (PJ), with obvious notation, as 



T^(M 2 (X),{MX)})^T^(M 2 (X),{fn(X)}) 

oo 

+ J2({MX)}-{M(X)})' 



s({M(x)}y 



(D5) 



Since, the difference x( x ) = M 2 (x) — m 2 is 0(A) [Q ffl] f° r a per turbative [loop] expansion, it is obvious that 
{MX)} — {/n(X)} = 0(A) [0(6)]. Substituting (D3) - (D5) into ( |D2| ), we obtain a new series. Truncating this series 
at the nth order, we obtain a perturbation series in theory II. It is obvious from the above argument that this scries 
is free from UV divergences. 



APPENDIX E: DERIVATIONS OF (8.53) - (8.55) 



Here we compute 

F(Xo,Y ;r) 



dp 



p sin(pr) 

p z — L0 Z 



sin 



(y^ 2- ^(Xv + F )) - sin (y/F=Z?(X Q - F )) 



\Jp 2 — UJ 2 



-2Y cos (y/p 2 -LU 2 (X Q - Y ) 



dp 



-2iY a 



p 2 — UJ 2 



( P 2 - uj 2 y>i 2 



p, a t— ± 

ppi 



p, a— ± 



^ per exp ppr + er\/p 2 — ^ 2 (Xo + tY ) J 
^ + a^/p 2 -LO 2 (X -Y )Yj 



In the above equations, Xp > Yq and u) should be understood to be w + if) (r/ = + ). 
Now we pick out from (|E^), 



f(r, Z ) 



„ ^(^-^)3/ 2 eX1) 



f-i pr- \Jp 2 - uS 2 Z ^ 



We note that, for r > Z n 



dp- 



r exp i 



pr - 



\/p 2 -uj 2 Z c 



= 0, 



(El) 



(E2) 



(E3) 



(E4) 



T C - (p 2 -u 2 Y /2 

where C_ is the closed contour in a complex-p plane: 

C- = -> i? -> -iR -► (i? ^ oo). 

Here the contour segment (0 — ► R) is on the real axis, (R — > —iR) is on Re~ l6 (6 = — > 7r/2), and (— i_R — > 0) is on 
the imaginary axis. The contribution from the contour segment (R — > — ii?) to (E4) vanishes in the limit i? — > oo. 
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The contribution from (—i R —* 0) is pure imaginary. Then, f(r, Zq) in (E3) is pure imaginary. Since (|E2| ) is pure real, 
when ( |E3| ) is used in (E^), (E3) should be ca nce led by other contribution(s), so that we can regard (|E3D as vanishing. 

In the complex-p plane, the integrand of (p§) has a branch point at p — u> + irj. We take the branch cut to be a 
straight line [uj + irj, — oo + irj). Then, for r < Zq, we have 



dp- 



T C+ 

where the closed contour C+ is defined as 



exp 



( — i pr - x/p 2 - u 2 Z Q ^ 



0. 



(E5) 



C+ = R iR -> i(r) + 0+) -> i(r? + 0+) + w - a 

-> i(»7 - 0+) + w - a -> i(r? - 0+) -> (iJ -► oo). 

Here, (i? — » iR) is on i?e l9 (0 = — * tt/2) and + + ) + cj — a — > i(r/ — + ) + a; — a) is on a circle of radius a, 
whose center is at p = ui + irj. From (fEq), we obtain, in the limit r\ — > + , 



f(r,Z ) = -2i dp 



pe 



,, r , - cosh( vV - P 2 Z ) + id — e-^ r - ttZqb 

q — p i yi A v ijjq. 



Here, as above, a pure- imaginary contribution has been dropped. Taking the limit a — > + , we obtain 



/(r,Z ) = -t- / ^ 



Wo s (l-e 2 ) 3/2 
2 



COS 



+i— e 



Another type of integral in (E2) may be computed similarly. Up to a pure-imaginary contribution, we obtain 

-i[pr-^/p 2 ~ui 2 Zq] 



g(r, Z ) = / dp 



9 2 



fl(Z - r) 



di ^— — smh (yi-e 2 ^o 



With these preliminaries, we can now compute (E2): 

T(X Q , Yq- r) = 9(X + Yq-t) 0(r - (X - Y )) - { ^(X + Y ) cos(wr) - sin(^r) 



o d? (l-^ 2 ) 3/2 



sin(wr£) cosh (y/ 1 — £ 2 w(Jfo + Yq)J — sin(wr) 



+0(X O - Y a - r) 



1 J f 1 £ sin(a;r£) 



s (w 2 ) 3/2 



T = ± 

-^u /' rf ^^^ sinh (VT~ecj(X Y ) 



Tcosh (x/l-eu{X + tY 







i-e 2 



It is straightforward to obtain the approximate form for J-(Xq, Yq; r) in several limits for the arguments: 
1. Xq-Yq <r < Xq + Yq; ujr < uj(X + Y ) « 1: 



F(X 07 Y ;r) ~ - 



w 2 (x + y ) s 



2. X - Yq < r < X + Y ; 1 << u(X Q + Y ) and tur « xfuj(X + Y ) 

/27T r 



F(X ,Y ;r) 



4 (o;(Z + Y )f/ 2 

7T 

+ -(^o + F )coswr, 



„o((X +Yb) 



h(u(X Q + Y (i ))-\-^— 



(E6) 
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where 



3. r<X Q -Y Q ;uj(X Q + Y Q )«l: 



Kx) = 1 + \~ + 
8 x 



HXo,Y ;r) ~ -(luY q ) 2 r. 



4.r<X -Y ; uj(X - Y ) « 1 « u(X + Y ): 



F(X ,Y ;r) 



2tt 



4 (U>(X + y ))3/a 



^(Xo+Yo) 



2 x + r 
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